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Abstract In this lectures* I discuss the electronic liquid crystal (ELC) phases 
in correlated electronic systems, what these phases are and in what context 
they arise. I will go over the strongest experimental evidence for these phases 
in a variety of systems: the two-dimensional electron gas (2DEG) in magnetic 
fields, the bilayer material Sr3Ru2 07 (also in magnetic fields), and a set of 
phenomena in the cupratc superconductors (and more recently in the pnictide 
materials) that can be most simply understood in terms of ELC phases. 
Finally we will go over the theory of these phases, focusing on effective field 
theory descriptions and some of the known mechanisms that may give rise to 
these phases in specific models. 



1 Electronic Liquid Crystal Phases 

Electronic liquid crystal phases [TT5] are states of correlated quantum elec- 
tronic systems that break spontaneously either rotational invariance or trans- 
lation invariance. Since most correlated electronic systems arise in a solid 
state environment the underlying crystal symmetry plays a role as it is the 
unbroken symmetry of the system. Thus in practice these phases break the 
point group symmetry of the underlying lattice, in addition of the possible 
breaking of the lattice translation symmetry. 

This point of view is commonplace in the classification of phases of clas- 
sical liquid crystals. [35] Classical liquid crystal systems are assemblies of a 
macroscopically large number of molecules with various shapes. The shapes 
of the individual molecules (the "nematogens" ) affects their mutual interac- 
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tions, as well as enhancing entropically-driven interactions ("steric forces") 
which, when combined, give rise to the dazzling phase diagrams of liquid 
crystals and the fascinating properties of their phases. [15] 

The physics of liquid crystals is normally regarded as part of "soft" con- 
densed matter physics, while the physics of correlated electrons is usually 
classified as part of "hard" condensed matter physics. The necessity to use 
both points of view clearly brings to the fore the underlying unity of Physics 
as a science. Thus, one may think of this area as "soft quantum matter" or 
quantum soft matter" depending on to which tribe you belong to. 

These lectures are organized as follows. In section [T]ELC phases and their 
symmetries are described. In section [2] I cover the main experimental evi- 
dence for these phases in 2DEGs, in Sr3Ru207 and in high temperature 
superconductors. In section [3] I present the theories of stripe phases, in sec- 
tion m the relation between electronic inhomogeneity and high temperature 
superconductivity is discussed, and sectioii5]is devoted to the theory of the 
pair density wave (the striped superconductor). Section |6] is devoted to the 
theories of nematic phases and a theory of nematic electronic order in the 
strong coupling regime is discussed in section [71 The stripe- nematic quantum 
phase transition is discussed in section [H 



1.1 Symmetries of Electronic Liquid Crystal Phases 
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Fig. 1 Cartoon of liquid crystal phases. 
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We will follow Rcf. |115] and classify the ELC phases of strongly eorre- 
lated electrons^ following the symmetry-based seheme used in classical liquid 
crystals [IS [23: 

1. Crystalline phases: phases that break all continuous translation symme- 
tries and rotational invariance. 

2. Smectic ("stripe") phases: phases break one translation symmetry and 
rotational invariance. 

3. Hematic and hexatic phases: uniform (liquid) phases that break rotational 
invariance. 

4. Isotropic: uniform and isotropic phases. 

A cartoon of the real space structure of these ordered phases is shown in 



Unlike classical liquid crystals, electronic systems carry charge and spin, 
and have strong quantum mechanical effects (particularly in the strong corre- 
lation regime) . This leads to a host of interesting possibilities of ordered states 
in which the liquid crystalline character of the spatial structure of these states 
becomes intertwined with the "internal" degrees of freedom of electronic sys- 
tems. These novel ordered phase will be the focus of these lectures. One of 
the aspects we will explore is the structure of their phase diagrams. Thus in 
addition of considering the thermal melting of these phases, we will also be 
interested in the quantum melting of these states and the associated quantum 
phase transitions, (see a sketch in FigEl) 

In this context, the crystalline phases are either insulating or "almost in- 
sulating", e.g. multiple charge density density waves (CDW) ordered states 
either commensurate or sliding (incommensurate). However, these phase may 
also be superconducting either by coexistence or, more interestingly, by mod- 
ulating the superconducting states themselves. Similarly, electron nematics 
are anisotropic metallic or superconducting states, while the isotropic phase 
are also either metallic or superconducting. As we will see these phases dis- 
play a set of rather striking and unusual behaviors, some of which have been 
observed in recent experiments. 



1.2 Order Parameters and their Symmetries 

The order parameters of ELC phases are well known, 151 1117) In the crys- 
talline phases, the order parameters are p Ki the expectation values of the 
density operators at the set of ordering wave vectors {K} the defines the 
crystal, [in] 



Fig[Il 




(1) 



^ You may call the ELC phases the anisotropic states of point particles! 
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Fig. 2 Schematic phase diagram of electronic liquid crystal phases. Here T is tem- 
perature and r denotes a tuning parameter the controls the strength of the quantum 
fluctuations. In practice it can represent doping, magnetic field, pressure or even 
material. The full dots are quantum critical points. 



where p(r) is the local charge density. Thus, under an uniform translation by 
Pk transforms as 

PK ^ PK ^ (2) 

Smcctic phases are unidirectional density waves and their order parameters 
are also expectation values pK but for only one wave vector K. For charged 
systems, p(r) is the charge density, and the order parameter pK is the charge 
density wave (CDW) order parameter. Since p{r) is real, pK = pIk' ^^^^ the 
density can be expanded as 

p(r)=po(r) + PK(r)e^^-'- + c.c. (3) 

where po{r) are the Fourier components close to zero wave vector, k = 0, and 
PK(r) are the Fourier components with wave vectors close to k = K. Hence, 
a density wave (a smectic) is represented by an complex order parameter 
field, in this case PK(r). This is how we will describe a CDW and a charge 
stripe (which from the point of view of symmetry breaking have the same 
description).^ 



^ On the other hand, in the case of a crystal phase, the expansion is 

p(r) = po (r) + E PK (r)e»^-- + c.c. (4) 

where F denotes the set of primitive lattice vectors of the crystal phase. [21] 
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Smcctic order is detected most easily in scattering experiments through 
the measurement of the static structure factor, usuaUy denoted by S'(k), 

Sik) ^ I ^SiKuj) (5) 

where SCk^uj) is the dynamical structure factor, i.e. the Fourier transform 
of the (in this case) density-density correlation function. The signature of 
smectic order is the existence of delta-function component of S'(k) at the 
ordering wave vector, k = K, with a prefactor that is equal to |(/3k)P-[1S] 

In the case of a spin density wave (a "spin stripe" ) the picture is the same 
except that the order parameter field is multi-component, SK(r), correspond- 
ing to different spin polarizations. Thus, the local spin density S(r) has the 
expansion 

S(r) = So(r) + SK(r)e'K---t-c.c. (6) 

where So(r) denotes the local (real) ferromagnetic order parameter and SK(r) 
is the (complex) SDW (or spin stripe) order parameter field, a complex vector 
in spin space. 

One of the questions we will want to address is the connection between 
these orders and superconductivity. The superconducting order parameter, 
a pair condensate, is the complex field ^(r). It is natural (and as we will 
see borne out by current experiments) to consider the case in which the 
superconducting order is also modulated, and admits an expansion of the 
form 

Air) = Z\o(r) + A^^ir) e^'^ + A^^{v) er'"^-- (7) 

where the uniform component Aq is the familiar BCS order parameter, and 
A-k{y) is the pair- density-wave (PDW) order parameter. [T?[ [Tl] closely re- 
lated to the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) order parameter [69l 
I124j (but without an external magnetic field). Since A{y) is complex, Aii_{v) ^ 
Z\_K(r)*, the PDW state has two complex order parameters.^ 

In contrast, nematic phases are translationally invariant but break rota- 
tional invariance. Their order parameters transform irreducibly under the 
rotation group for a continuous system, or under the point (or space) group 
of the lattice. Hence, the order parameters of a nematic phase (hexatic and 
their generalizations) are symmetric traceless tensors, that we will denote by 
Qij ■ [45] In 2D as most of the problems we will be interested in are 2D systems 
(or quasi2D systems), the order parameter takes the form (with i,j = x,y) 

which, alternatively, can be written in terms of a director N, 



^ I will not discuss the case of spiral order here. 
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xy — 



\N\ e*^ 



(9) 



Under a rotation by a fixed angle 6, N transforms as* 



N ^ N e' 



(10) 



Hence, it changes sign under a rotation by tt/2 and it is invariant under 
a rotation by tt (hence the name director, a headless vector). On the other 
hand, it is invariant under uniform translations by R. 

In practice we will have a great latitude when choosing a nematic order 
parameter as any symmetric traceless tensor in space coordinates will trans- 
form properly. In the case of a charged metallic system, a natural choice to 
describe a metallic nematic state is the traceless symmetric component of the 
resistivity (or conductivity) tensor. [551 [551 [51 [T ^ In 2D we will use 



where p^x and pyy are the longitudinal resistivities and Pxy = Pyx is the 
transverse (Hall) resistivity. A similar analysis can be done in terms of the 
dielectric tensor, which is useful in the context of light scattering experi- 
ments. On the other hand, when looking at the spin polarization properties 
of a system other measures of nematic order arc available. For instance, in 
a neutron scattering experiment, the anisotropy under a rotation TZ (say, by 
7r/2) of the the structure factor S(k) 



is a measure of the nematic order parameter Q. [1171 [90] . 

Other, more complex, yet quite interesting phases are possible. One should 
keep in mind that the nematic order parameter (as defined above) corresponds 
to a field that transforms under the lowest (angular momentum £ = 2) irre- 
ducible representation of the group rotations, compatible with inversion sym- 
metry. The nematic phase thus defined has d-wavc symmetry, the symmetry 
of a quadrupole. Higher symmetries arc also possible, e.g. hexatic {i — 6). 
However it is also possible to have states that break both rotational invariance 
and 2D inversion (mirror refiection), as in the ^ = 3 channel. Such states break 
(although mildly) time- reversal invariance. j 1941 r206| Other complex phases 
arise by combining the nematic order in real space with those of some internal 
symmetry, e.g. spin or orbital degeneracies. Thus one can consider nematic 
order parameters in the spin-triplet channels, which give rise to a host of 
(as yet undetected) phases with fascinating behaviors in the spin channel or 




Pxx Pyy 
Pxy 



Pxy 
Pyy Pxx 



(11) 



Q - S'(k) - S'(7^k) 



(12) 



* For a lattice system, rotational symmetries are those of the point (or space) group 
symmetry of the lattice, thus, nematic order parameters typically become Ising-like 
(on a square lattice) or three-state Potts on a triangular lattice (and so forth.) 
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under timc-rcversal, such as the dynamical generation of spin orbit couphng 
or the spontaneous breaking of time- reversal invariance. |220[ [2l9l 1194] 



1.3 Electronic Liquid Crystal Phases and Strong 
Correlation Physics 

One of the central problems in condensed matter physics is the understanding 
of doped Mott insulators. Most of the interesting systems in condensed mat- 
ter, notably high temperature superconductors, are doped Mott insulators. [3] 
A Mott insulator is a phase of an electronic system in which there is a gap 
in the single particle spectrum due to the effects of electronic correlations 
and not to features of the band structure. Thus, Mott insulators have an odd 
number of electrons in the unit cell. For a system like this, band theory would 
predict that such systems must be metallic, not insulating, and be described 
by the Landau theory of the Fermi liquid. Electronic systems that become 
insulating due to the effects of strong correlation are states of matter with 
non-trivial correlations. 

Most known Mott insulating states are ordered phases, associated with 
the spontaneous breaking of some global symmetry of the electronic system, 
and have a clearly defined order parameter. Typically the Mott state is an 
antiferromagnetic state (or generalizations thereof). However there has been 
a sustained interest possible in non-magnetic Mott phases, e.g. dimerized, 
various sorts of conjectured spin liquids, etc., some of which do not admit an 
order parameter description (as in the case of the topological phases.) 

We will not concern ourselves on these questions here. What will matter 
to us is that doping this insulator by holes disrupts the correlations that 
define the insulating state. Consequently doped holes are more costly (ener- 
getically) if they are apart than if they are together. The net effect is that 
the disruption of the correlations of the Mott state results in an effective 
strong attractive interaction between the doped holes. This effect was early 
on mistaken for aa sign of pairing in models of high temperature supercon- 
ductors (such as the Hubbard and t — J models). Further analysis revealed 
that this effective attraction meant instead the existence of a generic insta- 
bility of strongly correlated systems to phase separation.\bb\ This feature of 
strong correlation has been amply documented in numerical simulations (see, 
for instance, Ref.pT]). 

Due to the inherent tendency to phase separation of Hubbard-type models 
(and its descendants), the insulating nature of a Mott insulator cannot be 
ignored and, in particular, its inability to screen the longer range Coulomb in- 
teractions. Thus, quite generally, one can expect that the combined effects of 
the kinetic energy of the doped holes and the repulsive Coulomb interactions 
should in effect frustrate the tendency to phase separation of short-ranged 
models of strong correlation. [54] 
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Fig. 3 Schematic phase diagram of the cuprate superconductors. The full lines are 
the phase boundaries for the antiferromagnetic and superconducting phases. The bro- 
ken line is the phase diagram for a system with static stripe order and a pronounced 
1/8 anomaly. The dotted line marks the crossover between the bad metal and pseu- 
dogap regimes. 

The existence of strong short range attractive forces and long range repul- 
sion is a recipe for the formation of phases with complex spatial structure. 
As noted above, this is what happens in classical liquid crystals. It is also the 
general mechanism giving rise to generally inhomogeneous phases in classi- 
cal complex fluids such as ferrofluids and heteropolvmers. [189| as well as in 
astrophysical problems such as the crusts of neutron stars. [T40] 

The point of view that we take in these lectures is that the behavior ob- 
served in the underdoped regime of high temperature superconductors and in 
other strongly correlated systems is due to the strong tendency that these sys- 
tems have to form generally inhomogeneous and anisotropic phases, "stripes" . 
In the following lectures we will go over the experimental evidence for these 
phases and for their theoretical underpinning.^ 



2 Experimental Evidence in Strongly Correlated 
Systems 

During the past decade or so experimental evidence has been mounting of 
the existence of electronic liquid crystal phases in a variety of strongly cor- 
related (as well as not as strongly correlated) electronic systems. Wc will be 
particularly interested in the experiments in the copper oxide high tempera- 



^ Ref. [64| is a recent, complementary, review of the phenomenology of nematic phases 
in strongly correlated systems. 
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Fig. 4 a) 2DEG in a magnetic field, b) Landau levels. 



ture superconductors, in the ruthcnate materials (notably Sr3Ru207), and in 
two-dimensional electron gases (2DEG) in large magnetic fields. However, as 
we will discuss below, these concepts are also relevant to more conventional 
CDW systems. 



2.1 Nematic Phases in the 2DEG in High Magnetic 
Fields 

To this date, the best documented electron nematic state is the anisotropic 
compressible state observed in 2DEGs in large magnetic fields near the middle 
of a Landau level, with Landau index N > 2[T3Sl[I31[5ni[IS5] (FigE]). In uhra 
high mobility samples of a 2DEG in AlAs-GaAs heterostructures, transport 
experiments in the second Landau level (and above) near the center of the 
Landau level show a pronounced anisotropy of the longitudinal resistance 
rising sharply below T ~ 80 mK, with an anisotropy that increases by orders 
of magnitude as the temperature is lowered. This effect is only seen in ultra- 
clean samples, with nominal mean free paths of about 0.5 mm (!) and nominal 
mobilities of 10 - 30 x 10^.^ 

A nematic order parameter can be constructed phenomenologically from 
the measured resistivity tensor, by taking the symmetric traceless piece of it. 
This was done in Ref.[BS] where a fit of the data of Lilly et al |1351 1134] was 



^ The anisotropy is strongly suppressed by disorder 
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Fig. 5 Spontaneous magneto-transport anisotropy in tlie 2DEG: a) peaks in p^a; 
developing at low T in liigli LLs (dotted line: T=100mK: thick line: 65mK; thin line: 
25mK). Inset: temperature dependence of peak height at u = 9/2 (closed circles), 
11/2 (open circles), 13/2 (closed triangles) and 15/2 (open triangles), b) Anisotropy 
of p^^ at T=25mK. (Lilly et alp^l. 



shown to be consistent with a classical 2D XY model (in a weak symmetry 
breaking field). A 2D XY symmetry is expected for a planar nematic order 





0.8 






+^ 


0.6 


Cl 






0.4 


Cl 






0.2 








■ Experiinental dala 
• h/J=0.05 J=73.5 mK 




100 200 

T(mK) 

(a) Left 



300 



l.^ 1.4 1.5 

B<T) 



a) 








V = 9/2 


1 




1 — R„ 






















1 




1 


n 



50 100 

T(Q!K) 




6 3 _ 3 S 
B„(T)[im] B,<T) [UOJ 



(b) Right 



Fig. 6 Left: Nematic order in the 2DEG; fit of the resistance anisotropy to a 2D 
XY model Monte carlo simulation. [65]. Right: a) Longitudinal resistance anisotropy 
around i/ = 9/2 at T = 25 mK. Solid trace: Rxx', average current flow along [110]. b) 
Temperature dependence of resistances at u = 9/2. c) R^x and Ryy at u = 9/2 at 
T = 25 mK vs in-plane magnetic field along [110] and [110]. (Cooper et al[38]'). 



provided the weak lattice symmetry breaking is ignored. Presumably lattice 
anisotropy is responsible for the saturation shown at low temperatures in 
FigEKleft panel). 
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These experiments were originally interpreted as evidence for a quan- 
tum Hall smectic (stripe) phase [1^01 11551 11441 [TT] . However, further 
experiments (|371 [351 ISS] (FiglH Right panel) did not show any evidence 
of pinning of this putative unidirectional CDW as the I ~ V curves were 
found to be strictly linear at low bias. In addition broadband noise in 
the current, a characteristic of CDW systems. In contrast, extremely sharp 
threshold electric fields and broadband noise in transport was observed in a 
nearby reentrant integer quantum Hall phase, suggesting a crystallized elec- 
tronic state. These facts, together with a detailed analysis of the experi- 
mental data, suggested that the compressible state is in an electron nematic 
phase [63l l65l 12141 11771 148] . which is better understood as a quantum melted 
stripe phase.'' An alternative picture, a nematic phase accessed by a Pomer- 
anchuk instability from a "composite fermion" Fermi liquid is conceivable 
but hard to justify microscopicallv. [1631 H5] 



2.2 The Nematic Phase of Sr^Ru^O^ 



1.4 




Field (tesia) 



Fig. 7 Phase diagram of Sr3Ru207 in the temperature-magnetic field plane, (from 
Grigera et al|73p. 



''' The 2DEG in a strong magnetic field is inherently a strongly correlated system as 
the interaction is always much bigger than the (vanishing) kinetic energy. 



12 



Eduardo Fradkin 



Recent magneto-transport experiments in the quasi-two-dimensional bi- 
layer ruthenate Sr3Ru207 by the St. Andrews group[T^ have given strong 
evidence of a strong temperature-dependent in-plane transport anisotropy in 
strongly correlated materials at low temperatures T < 800 mK and for a win- 
dow of perpendicular magnetic fields around 7.5 Tesla.(Se Figl?!) Sr3Ru207 
is a quasi-2D bilayer material known to have a metamagnetic transition as a 
function of applied perpendicular magnetic field and temperature.® As in the 
case of the 2DEG of the previous section, the transport anisotropy appears at 
very low temperatures and only in the cleanest samples. The observed trans- 
port anisotropy has a strong temperature and field dependence (although not 
as pronounced as in the case of the 2DEG) is shown in Fig|8] These exper- 
iments provide strong evidence that the system is in an electronic nematic 
phase in that range of magnetic fields |19[ I62| . The electronic nematic phase 
appears to have preempted a metamagnetic QCP in the same range of mag- 
netic fields [711 1152[|166[ [7^ . This suggests that proximity to phase-separation 
may be a possible microscopic mechanism to trigger such quantum phase 
transitions, consistent with recent ideas on the role of Coulomb-frustrated 
phase separation in 2DEGs piI51 1141| . 




Fig. 8 Left: paa and pi,;, of the in-plane magnetoresistivity tensor of a high-purity 
single crystal of Sr3Ru207. (A) For an applied c-axis field, paa (black) and pj,;, 
(red). (B) Tilted field (13° from c). Right: The temperature dependence difference 
of paa and Pi,;, for fields applied at 6 = 72° such that the in-plane field component 
lies along a. (Upper inset) Temperature dependence of paa (black) and p(,{, (red) for 
Ij,qH = 7.4 T applied in the direction specified above. (Lower inset) The temperature 
dependence of the difference between the two magnetoresistivities shown in the upper 
inset, normalized by their sum. (from Borzi et al|19|). 



® Contrary to the case of the 2DEG in AlAs-GaAs heterostructures and quantum 
wells, the magnetic fields applied to Sr3Ru207 are too weak to produce Landau 
quantization. 
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2.3 Stripe Phases and Nematic Phases in the Cuprates 



In addition to high temperature superconductivity, the copper oxides mate- 
rials display a strong tendency to have charge-ordered states, such as stripes. 
The relation between charge ordered states [TT4] . as well as other proposed 
ordered states[30 ll206j . and the mechanisni(s) of high temperature supercon- 
ductivity is a subject of intense current research. It is not, however, the main 
focus of these lectures. 
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Fig. 9 Left: Static spin and charge stripe order in La2-a;Ba2,Cu04 in neutron scat- 
tering (From Fujita et al|68p. Right: resonant X-ray scattering (From Abbamonte et. 



Stripe phases have been extensively investigated in high temperature su- 
perconductors and detailed and recent reviews are available on this subject [11 71 
I202j . Stripe phases in high temperature superconductors have unidirectional 
order in both spin and charge (although not always) which are typically in- 
commensurate. In general the detected stripe order (by low energy inelastic 
neutron scattering) in La2-a;Sra;Cu04, La2-a;Ba2;Cu04 and YBa2Cu306+a; 
(see Rcfs. |117] and |202| and references therein) is not static but "fluctu- 
ating". As emphasized in Ref. [117j . "fluctuating order" means that there is 
no true long range unidirectional order. Instead, the system is in a (quan- 
tum) disordered phase, very close to a quantum phase transition to such an 
ordered phase, with very low energy fluctuations that reveal the character 
of the proximate ordered state. On the other hand, in La2-2:Baa;Cu04 near 
j; = 1/8 (and in Lai.6-a;Ndo.4Sra;Cu04 also near x = 1/8 where they were dis- 
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Fig. 10 Nematic order in underdoped YBa2Cu306+a; [y = 6.45). (a to c) Intensity 
maps of the spin-excitation spectrum at 3, 7, and 50 meV, respectively. The a* and 
6* directions are indicated in (a), (d) Colormap of the intensity at 3 meV, as it would 
be observed in a crystal consisting of two perpendicular twin domains with equal 
population, (e and f) Scans along a* and 6* through Q^^. (from Hinkov et al|90|). 

covered first |200| ). the order detected by elastic neutron scattering |201| . and 
resonant X-ray scattering in La2-xBaxCu04 jl] also near x = 1/8, becomes 
true long range static order. (See FiglS]) 

In the case of La2-xSra;Cu04, away from x = 1/8, and particularly 
on the more underdoped side, the in-plane resistivity has a considerable 
temperature-dependent anisotropy, which has been interpreted as an indi- 
cation of electronic nematic order. [5] The same series of experiments also 
showed that very underdoped YBa2Cu306-i-a; is an electron nematic as well. 

The most striking evidence for electronic nematic order in high tem- 
perature superconductors are the recent neutron scattering experiments 
in YBa2Cu306+:r at y = 6.45[88j. (See Figs. [IO][ni)-) In particular, the 
temperature-dependent anisotropy of the inelastic neutron scattering in 
YBa2Cu306+i^ shows that there is a critical temperature for nematic order 
(with Tc ^ 150 A") where the inelastic neutron peaks also become incom- 
mensurate. Similar effects were reported by the same group [89] at higher 
doping levels {y 6.6) who observed that the nematic signal was decreasing 
in strength suggesting the existence of a nematic-isotropic quantum phase 
transition closer to optimal doping. Fluctuating stripe order in underdoped 
YBa2Cu306+z: has been detected earlier on in inelastic neutron scattering 
experiments |156i 1193] which, in hindsight, can be reinterpreted as evidence 
for nematic order. However, as doping increases past a j/ ~ 6.6 a spin gap ap- 
pears and magnetic scattering is strongly suppressed at low energies (in the 
absence of magnetic fields) making inelastic neutron scattering experiments 
less effective in this regime. 

In a particularly series of interesting experiments, the Nernst coefficient 
was measured in YBa2Cu306+x ranging from the very underdoped regime. 
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Fig. 11 a) Incommensurability 6 (red symbols), half- width-at-half- maximum of the 
incommensurate peaks along a* ^ , black symbols) and along b* {^^ ^ , open blue 
symbols) in reciprocal lattice units, (from Hinkov et al[90|'). b) Static stripe order 
induced by an external magnetic field in YBa2Cu3 06+x at y = 6.45. 



where inelastic neutron scattering detects nematic order, to a slightly over- 
doped regime. [44] 

The Nernst coefficient is defined as follows. Let je and jg be the charge 
and heat currents established in a 2D sample by an electric field E and a 
temperature gradient V^- 

where a, a and k are 2x2 tensors for the conductivity, the thermoelectric 
conductivity and the thermal conductivity respectively. The Nernst coeffi- 
cient, also a 2 X 2 tensor 6 is measured (see Ref. [132j ) say by applying a 
temperature gradient in the x direction and measuring the voltage along the 
y direction: 

E = -9\jT (14) 
Since no current flows through the system, the Nernst tensor is 

e = -a-^a (15) 

These experiments revealed that the Nernst (tensor) coefficient has an 
anisotropic component whose onset coincides (within the error bars) with 
the conventionally defined value of the pseudogap temperature T*, and es- 
sentially tracks its evolution as a function of doping. Thus, it appears that. 
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at least in YBa2Cu306+a;, the pseudogap is a regime with nematic order (see 
Figlll)9 
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Fig. 12 The pseudogap as a regime with of charge nematic order in YBa2Cu306+a;: 
measured from the Nernst effect, the onset coincides with tlie T* of the pseudogap. 
(From Ref.gl] and Rei.^7\.) 

Inelastic neutron scattering experiments have found nematic order also in 
La2_a;Sra;Cu04 materials where fluctuating stripes where in fact first discovered |200j . 
Matsuda et al }145j have found in undcrdopcd La2-a:Sr2:Cu04 (x — 0.05), a 
material that was known to have "fluctuating diagonal stripes", evidence 
for nematic order similar to what Hinkov et al[HBl found in underdoped 
YBa2Cu306+2^. Earlier experiments in La2-a;Sra;Cu04 in moderate magnetic 
fields had also shown that a spin stripe state became static over some critical 
value of the field. [7] These experiments strongly suggest that the experiments 
that had previously identified the high temperature superconductors as hav- 
ing "fiuctuating stripe order" (both inside and outside the superconducting 
phase) were most likely detecting an electronic nematic phase, quite close to 
a state with long range stripe (smectic) order. In all cases the background 
anisotropy (due to the orthorhombic distortion of the crystal structure) acts 
as a symmetry breaking field that couples linearly to the nematic order, 
thus rounding the putative thermodynamic transition to a state with sponta- 
neously broken point group symmetry. These effects are much more apparent 
at low doping where the crystal orthorhombicity is significantly weaker. 

In La2-a;Baa;Cu04 at a; = 1/8 there is strong evidence for a complex stripe 
ordered state that intertwines charge, spin and superconducting order [T33l ITS ] 
(shown in FigHS]) In fact La2-a;Ba2:Cu04 at x = 1/8 appears to have some 
rather fascinating properties. As summarized in FiglHl La2-a;Baa;Cu04 at 
a; = 1/8 has a very low critical superconducting Tc ^ 4K (where the Meiss- 



^ The same group had shown earlier than the Nernst coefficient is a sensitive indicator 
of the onset of stripe charge order in Lai R-^Ndn ASr^CuO^|40| 
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Fig. 13 

data, b) 
al[T33]). 



a) Dynamical layer decoupling in La2-a;Baa;Cu04 at a: = 1/8 from transport 
Kosterlitz-Thouless transition in La2_a;Baa;Cu04 at a; = 1/8 (from Q. Li et 



54 K .. 



T ■ = 42 K 

spin 



T"=35K 



T„M6K 



T, = 4K 
OK 




max {p j-8xlO-mi2 cm 



P. ~ p _ lOmfl 



p - p - 

ab c 

Meissner State 



Fig. 14 Summary of the behavior of the stripe-ordered superconductor 
La2-a;Baa;Cu04 near x = 1/8: Tco is the charge ordering temperature, Tgpin the 
spin ordering temperature, T** marks the beginning of layer decoupling behavior, 
Tkt is the 2D superconducting temperature ("KT"), T^d is the 3D resistive transi- 
tion, and Tc is the 3D Meissner transition, (from Ref.|15|) 
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ner state sets in). However it is known from angle-resolved photoemission 
(ARPES) experiments that the anti-nodal gap (which roughly gives the pair- 
ing scale) is actually largest at a; = 1/8 |204| (or unsuppressed by the 1/8 
anomaly according to Ref. [S^.) Static charge stripe order sets in at 54K 
(where there is a structural transition from the LTO to the LTT lattice struc- 
ture), but static spin stripe order only exists below 42 K. As soon as static 
spin order sets in, the in-plane resistivity begins to decrease very rapidly with 
decreasing temperature, while the c-axis resistivity increases. (See FigHSl Left 
panel). Below 35 K strong 2D superconducting fluctuations are observed and 
at 16 K the in-plane resistivity vanishes at what appears to be a Kosterlitz- 
Thouless transition, (shown in FigflBl right panel.) However, the full 3D 
resistive transition is only reached at 10 K (where Pc — >■ 0) although the 
Meissner state is only established below 4K! This dazzling set of phenom- 
ena shows clearly that spin, charge and superconducting order are forming 
a novel sort of intertwined state, rather than compete. We have conjectured 
that a pair density wave (PDW) is stabilized in this intermediate temperature 
regime [lll[T4j [16]. Similar phenomenology, i.e. a dynamical layer decoupling, 
has been seen in La2-a;Sra;Cu04 at moderate fields where the stripe order is 
static. [187) We will return below on how a novel state, the pair-density wave, 
explains these phenomena. 




(a) (b) 

Fig. 15 Short range stripe order in Bi2Sr2CaCu2 08+A- (Dy-Bi2212) as seen in STM 
experiments, (a) R maps taken at 150 mV; (b)v^i? shows the local nematic order, 
(from Kohsaka et al |118p . 

An important caveat to the analysis we presented here is that in doped sys- 
tems there is always quenched disorder, which has different degrees of short 
range "organization" in different high temperature superconductors. Since 
disorder also couples linearly to the charge order parameters it ultimately 
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also rounds the transitions and render the system to a glassy state (as noted 
already in Refs. jll51lll7| ). Such effects are evident in scanning tunneling mi- 
croscopy (STM) experiments in Bi2Sr2CaCu2 08+5 which revealed that the 
high energy (local) behavior of the high temperature superconductors has 
charge order and it is glassy [97l 11171 [83l 11181 1208] . This is most remarkable 
as the STM data on Bi2Sr2CaCu208+(5 at low bias shows quasiparticle prop- 
agation in the superconducting state (but not above Tc). Yet, at high bias 
(i.e. high energies) there are no propagating "quasiparticles" but, instead, 
provides a vivid image of electronic inhomogeneity with short range charge 
order. This behavior is contrary to what is commonly the case in conven- 
tional superconductors where STM at high energies shows Fermi-liquid like 
electronic quasiparticles. Similarly, the high energy spectrum of ARPES has 
never resembled that of a conventional metal. 

Finally, we note that in the recently discovered iron pnictides based family 
of high temperature superconductors, such as La (Oi_a;F2:)FeAs and and 
CafFei-^Co-r)9As9 |104[ I158j . a unidirectional spin-density- wave has been 
found. It has been suggested [60] that the undoped system LaOFeAs and 
CaFe2As2 may have a high-temperature nematic phase and that quantum 
phase transitions also occur as a function of fluorine doping [22T| [33 ] . This 
suggests that many of the ideas and results that we present here may be 
relevant to these still poorly understood materials. 

The existence of stripe-ordered phases is well established in other complex 
oxide materials, particularly the manganites and the nickelates. In general, 
these materials tend to be "less quantum mechanical" than the cuprates in 
that they are typically insulating (although with interesting magnetic proper- 
ties) and the observed charge-ordered phases are very robust. These materials 
typically have larger electron-phonon interactions and electronic correlation 
are comparatively less dominant in their physics. For this reason they tend to 
be "more classical" and less prone to quantum phase transitions. However, 
at least at the classical level, many of the issues we discussed above, such 
as the role of phase separation and Coulomb interactions, also play a key 
role[42]. The thermal melting of a stripe state to a nematic has been seen in 
the manganite material Bi^Ca^MnO':! [T84] . 



2.4 Conventional CDW materials 

CDWs have been extensively studied since the mid-seventies and there are 
extensive reviews on their properties |75l I76| . From the symmetry point of 
view there is no difference between a CDW and a stripe (or electron smec- 
tic). CDW states are usually observed in systems which are not particu- 
larly strongly correlated, such as the quasi-one-dimensional and quasi-two- 

A recent analysis of this data by Lawler and coworkers has revealed the existence of 
nematic order over much longer length scales than the broken positional order | 129| . 
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dimensional dichalcogenides, and the more recently studied tritellurides. 
These CDW states are reasonably well described as Fermi liquids (FL) which 
undergo a CDW transition, commensurate or incommensurate, triggered by 
a nesting condition of the FS |146[ [147] . As a result, a part or all of the FS is 
gapped in which case the CDW may or may not retain metallic properties. 
Instead, in a strongly correlated stripe state, which has the same symmetry 
breaking pattern, at high energy has Luttinger liquid behavior |1 151 153 } I25j. 

What will interest us here is that conventional quasi- 2D dichalcogenides, 
the also quasi-2D tritellurides and other similar CDW systems can quantum 
melt as a function of pressure in TiSe2 |192j , or by chemical intercalation as 
in Cu^TiSea [157113 and Nb^^TaSagS]. Thus, CDW phases in chalcogenides 
can serve as a weak-coupling version of the problem of quantum melting 
of a quantum smectic. Interestingly, there is strong experimental evidence 
that both TiSe2[Ill] and NbT.TaS9[15] do not melt directly to an isotropic 
Fermi fluid but go instead through an intermediate, possibly hexatic, phase. 
Whether or not the intermediate phases are anisotropic is not known as no 
transport data is yet available in the relevant regime. 

The case of the CDWs in tritellurides is more directly relevant to the 
theory we will present here. Tritellurides are quasi-2D materials which for a 
broad range of temperatures exhibit a unidirectional CDW {i.e. an electronic 
smectic phase) and whose anisotropic behavior appears to be primarily of 
electronic origin [22j 11251 11851 l59] . However, the quantum melting of this 
phase has not been observed yet. Theoretical studies have also suggested 
that it may be possible to have a quantum phase transition to a state with 
more than one CDW in these materials [257]. 



3 Theories of Stripe Phases 

3.1 Stripe Phases in Microscopic Models 

Of all the electronic liquid crystal phases, stripe states have been studied 
most. There are in fact a number of excellent reviews on this topic [25l 11171 
I210j covering both the phenomenology and microscopic mechanisms. I will 
only give a brief summary of important results and refer to the literature for 
details. 

As we noted in section II. 2[ stripe and CDW (and SDW) phases have 
the same order parameter as they correspond to the same broken symmetry 



Cua;TiSe2 is known to become superconducting j 1 57| . The temperature-pressure 
phase diagram of TiSe2 exhibits a superconducting dome enclosing the quantum 
critical point at whicli tlie CDW state melts. p"22] 
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state, and therefore the same order parameter pK (and Sq). ^'^ There is 
however a conceptual difference. CDW and SDW are normaUy regarded as 
weak coupUng instabiUties of a Fermi hquid (or free fermion state) typicahy 
triggered by a nesting condition satisfied by the ordering wave vector. pL 461 
I75j In this context, the quasiparticle spectrum is modified by the partial 
opening of gaps and a change in the topology of the original Fermi surface 
(or, equivalently, by the formation of "pockets"). Because of this inherently 
weak coupling physics, the ordering wave vector is rigidly tied to the Fermi 
wave vector kp. 

In one- dimensional systems, non-linearities lead to a more complex 
form of density wave order, a lattice of solitons, known in this context as 
discommensurations |147l I21| whose ordering wavevector is no longer neces- 
sarily tied to kp. A stripe state is essentially a two-dimensional generally 
incommensurate ordered state which is an analog of this strong coupling 
one-dimensional lattice of discommensurations . j 1 1 Oj Thus, in this picture, 
the spin stripe seen in neutron scattering [200] is pictured as regions of an- 
tiferromagnetic (commensurate) order separated by anti-phase domain walls 
(the discommensurations) where the majority of the doped holes reside. This 
picture is quite hard to achieve by any weak coupling approximation such as 
Hartree-Fock. 

Stripe phases were first found in Hartree-Fock studies of Hubbard and 
t — J models in two dimensions. |228[ I188[ 1168] In this picture stripe phases 
are unidirectional charge density waves (CDW) with or without an associated 
spin-density-wave (SDW) order. As such they are characterized by a CDW 
and/or SDW order parameters, pK and Sq respectively.^^ A Hartree-Fock 
theory of stripe phases was also developed in the context of the 2DEG in 
large magnetic fields [155] 11201 1144] to describe the observed and very large 
transport anisotropy we discussed above. 

As it is usually the case, a serious limitation of the Hartree-Fock approach 
is that it is inherently reliable only at weak coupling, and hence away from 
the regime of strong correlation of main interest. In particular, all Hartree- 
Fock treatments the stripe ground state typically produce an "empty stripe 
state", an insulating crystal and therefore not a metallic phase. Thus, in 
this approach a conducting (metallic) stripe phase can only arise from some 
sort of quantum melting of the insulting crystal and hence not describable 
in mean-field theory. The phenomenological significance of stripe phases was 
emphasized by several authors, particularly by Emery and Kivelson[54] 11111 
[56]. 

Mean field theory predicts that, at a fixed value of a electron (or hole) 
density (doping) , the generally incommensurate ordering wave vectors satisfy 
the relation K = 2Q. That this result should generally hold follows from a 
simple Landau-Ginzburg (LG) analysis of stripe phases (see, e.g. [1721 [T8] 

1^ In principle the order parameter of the stripe state may not be pure sinusoidal and 
will have higher harmonics of the fundamental order parameter. 

I will ignore here physically correct but more complex orders such as helical phases. 
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where it is easy to see that a trihncar term of the form p^^Q ' (^-nd its 
complex conjugate) is generahy aUowed in the LG free energy. In ar ordered 
state of this type the antifcrromagnctic spin order is "deformed" by anti- 
phase domain waUs whose periodic pattern coincides with that of the charge 
order, as suggested by the observed magnetic structure factor of the stripe 
state first discovered in the cuprate Lai s-rNdn ztSr^CuOA. [200] 

This pattern of CDW and SDW orders has suggested the popular cartoon 
of stripe phases as antiferromagnetic regions separated by narrow "rivers of 
charge" at antiphase domain walls. The picture of the stripe phase as an array 
of rivers suggests a description of stripe phases as a quasi-one-dimensional 
system. As we will see in the next subsection, this picture turned out to be 
quite useful for the construction of a strong coupling theory of the physics of 
the stripe phase. On the other hand, it should not be taken literally in the 
sense that these rivers always have a finite width which docs not have to be 
small compared with the stripe period and in many cases they may well be 
of similar magnitude. Thus, one may regard this phase as being described by 
narrow ID regions with significant transversal quantum fluctuations in shape 
(as it was presented in Ref. |115] ) or, equivalently, as quasi-lD regions with a 
significant transversal width. 

An alternative picture of the stripe phases can be gleaned from the t — J 
model, the strong coupling limit of the Hubbard model. Since in the result- 
ing effective model there is no small parameter, the only (known) way to 
treat it is to extend to SU{2) symmetry of the Hubbard (and Heisenberg) 
model to either SU{N) or Sp{N) and to use the large N expansion to study 
its properties. [MH [^TTl [^T^ [^TU] In this (large N) limit the undopcd anti- 
ferromagnet typically has a dimerized ground state, a periodic (crystalline) 
pattern of valence bond spin singlets. Since all degrees of freedom are bound 
into essentially local singlets this state is a quantum paramagnet. However, 
it is also "striped" in the sense that the valence bond crystalline state breaks 
at least the point group symmetry C4 of the square lattice as well as trans- 
lation invariance: it is usually a period 2 columnar state. In the doped 
system the valence bond crystal typically becomes a non-magnetic incom- 
mensurate insulating system. Mean- field analyses of these models [2T2] also 
suggest the existence of superconducting states, some of which are "striped" . 
Similar results are suggested by variational wave functions based of the RVB 
state. [131 |571 1224] We should note, however, that mean- field states are no 
longer controlled by a small parameter, such as and hence it is unclear 
how reliable they may by at the physically relevant case N ^ 2. For a detailed 
(and up-to-date) review of this approach see Ref. [210j . 

There are also extensive numerical studies of stripe phases in Hubbard type 
models. The best numerical data to date is the density matrix renormalization 
group (DMRG) work of White and Scalapino (and their collaborators) on 



This state is a close relative of the resonating valence bond (RVB) state originally 
proposed as a model system for a high superconducting stated] Ilt3 |. a (non- 
resonating) VB state. 
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Hubbard and i — J ladders of various widths (up to 5) and varying particle 
densities |215[|216l 12171 f^T5] and by Jackclmann et al in fairly wide ladders (up 
to 7).|77| An excellent summary and discussion on the results from various 
numerical results (as well as other insights) can be found in Ref.[lS]. The 
upshot of all the DMRG work is that there arc strong stripe correlations in 
Hubbard and t — J models which may well be the ground state. "'^^ 

Much of the work on microscopic mechanisms of stripe formation has been 
done in models with short range interactions such as the Hubbard and t — 
J models. As it is known. [1121 [551 HTj models of this type have a strong 
tendency to electronic phase separation. As we noted in the introduction, the 
physics of phase separation is essentially the disruption of the correlations 
of the the Mott (antiferromagnetic) state by the doped holes which leads 
to an effective attractive interaction among the charge carriers. When these 
effects overwhelm the stabilizing effects of the Fermi pressure {i.e. the fermion 
kinetic energy), phase separation follows. In more realistic models, however, 
longer range (and even Coulomb) interactions must be taken into account 
which tend to frustrate this tendency to phase separation [54], as well as a 
more complex electronic structure. |52j The structure of actual stripe phases 
in high Tc materials results from a combination of these effects. One of the 
(largely) unsolved questions is the relation between the stripe period and 
the filling fraction of each stripe at a given density. Most simple minded 
calculation yield simple commensurate filling fractions for each stripe leading 
to insulating states. At present there are no controlled calculations of these 
effects, although variational estimates have been published [142], except for 
results from DMRG studies on wide ladders. [77] 



3.2 Phases of Stripe States 

We will now discuss the strong coupling picture of the stripe phases. |53l 
I212[ [Til [H |25l 1114] . We will assume that a stripe phase exists with a fixed 
(generally incommensurate) wave vector K and a fixed filling fraction (or 
density) on each stripe. In this picture a stripe phase is equivalent to an 
array of ladders of certain width. In what follows we will assume that each 
stripe has a finite spin gap: a Luther-Emery liquid. [511 1143] 

3.2.1 Physics of the 2-leg ladder 

The assumption of the existence of a finite spin gap in ladders can be justified 
in several ways. In DMRG studies of Hubbard and t — J ladders in a rather 
broad density range, < a; < 0.3, it is found that the ground state has a finite 



A difficulty in interpreting the DMRG results lies in the boundary conditions that 
are used that tend to enhance inhomogeneous, stripe-like, phases. 
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spin gap. [162] Similar results were found analytically in the weak coupling 



Why there is a spin gap? There is actually a very simple argument for 
it.|57) In the non-interacting limit, U ~ V — 0, the two-leg ladder has two 
bands with two different Fermi wave vectors, ppi ^ pf2- Let us consider 
the effects of interactions in this weak coupling regime. The only allowed 
processes involve an even number of electrons. In this limit is is easy to see 
that the coupling of CDW fluctuations with Qi = 2pFi ^ Q2 — 2pF2 is 
suppressed due to the mismatch of their ordering wave vectors. In this case, 
scattering of electron pairs with zero center of mass momentum from one 
system to the other is a peturbatively (marginally) relevant interaction. The 
spin gap arises since the electrons can gain zero-point energy by delocalizing 
between the two bands. To do that, the electrons need to pair, which may 
cost some energy. When the energy gained by delocalizing between the two 
bands exceeds the energy cost of pairing, the system is driven to a spin-gap 
phase. 

This physics is borne out by detailed numerical (DMRG) calculations, even 
in systems with only repulsive interactions. Indeed, at a; = (the undoped 
ladder) the system is in a Mott insulating state, with a unique fully gapped 
ground state ( "CO^O" in the language of Ref. [5] ) ■ In the strong coupling limit 
(in which the rungs of the ladder are spin singlet valence bonds), U ^ t, the 
spin gap is large: As ~ J/2. |203j 

At low doping, < x < Xc ^ 0.3, the doped ladder is in a Luther-Emery 
liquid: there is no charge gap and large spin gap ("Cl^O"). In fact, in this 
regime the spin gap is found to decrease monotonically as doping increases. 
As 4, as x t, and vanishes at a critical value Xc'- As as a; Xc- 

The most straightforward way to describe this system is to use bosoniza- 
tion. Although in the ladder system has several bands of electrons that have 
charge and spin degrees of freedom, in the low energy regime the effective 
description is considerably simplified. Indeed, in this regime it is sufficient 
to consider only one effective bosonized charge field and one bosonized spin 
field. Since there is a spin gap, Ag ^ 0, the spin sector is massive. In contrast, 
the charge sector is only massive at x = 0, where there is a finite Mott gap 
Am. 

The effective Hamiltonian for the charge degrees of freedom in this 
(Luther-Emery) phase is 



where is the CDW phase field, and 9 is the SC phase field. They satisfy 
canonical commutation relations 



regime. [Sl fnTlfTMl [55] 



H = 



Jdy^ (5,0)Vi^(a.0)']+... 



(16) 



(17) 
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The parameters of this effective theory, the spin gap , the charge Luttinger 
parameter K, the charge velocity Vc, and the chemical potential n, have non- 
universal but smooth dependences on the doping x and on the parameters of 
the microscopic Hamiltonian, the hopping matrix elements t' /t and the Hub- 
bard interaction U /t. The ellipsis ... in the effective Hamiltonian represent 
cosine potentials responsible for the Mott gap Am in the undoped system 
(a; = 0). It can be shown that the spectrum in the low doping regime, a; — > 0, 
consists of gapless and spinless charge 2e fermionic solitons. 

The charge Luttinger parameter is found to approach A'— >l/2asa;— >0. 
As X increases, so does K reaching the value ~ 1 for x ^ 0.1. On the 
other hand if ~ 2 for x ^ Xc where the pin gap vanishes. The temperature 
dependence of the superconducting and CDW susceptibilities have the scaling 
behavior 

XSC ~ rp2-K 

(18) 
(19) 

Thus, both susceptibilities diverge Xcdw(T') ~^ oo and xsciT) — > oo for < 
X < Xc- However, for x < 0.1, the SC susceptibility is more divergent: xsc ^ 
XCDW- Hence, the doped ladder in the Luther-Emery regime is effectively a 
ID superconductor even for a system with nominally repulsive interactions 
(i.e. without "pairing"). 



3.2.2 The spin-gap stripe state 

Now consider a system with N stripes, each labeled by an integer a = 
1, . . . ,N. We will consider first the phase in which there is a spin gap. Here, 
the spin fluctuations are effectively frozen out at low energies. Nevertheless 
each stripe a has two degrees of freedom |115) : a transverse displacement field 
which describes the local dynamics of the configuration of each stripe, and 
the phase field 0a for the charge fluctuations on each stripe. The action of 
the generalized Luttinger liquid which describes the smcctic charged fluid of 
the stripe state is obtained by integrating out the local shape fluctuations as- 
sociated with the displacement fields. These fluctuations give rise to a flnite 
renormalization of the Luttinger parameter and velocity of each stripe. More 
importantly, the shape fluctuations, combined with the long- wavelength inter- 
stripe Coulomb interactions, induce inter-stripe density-density and current- 
current interactions, leading to an imaginary time Lagrangian density of the 
form 

/:smcctic = ^ .];{^)W^i<^~a')j-^{x). (20) 
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These operators arc marginal, i.e. have scahng dimension 2, and preserve the 
smectic symmetry (j)a 4>a + eta (where aa is constant on each stripe) of the 
decoupled Luttingcr fluids. Whenever this symmetry is exact, the charge- 
density- wave order parameters of the individual stripes do not lock with each 
other, and the charge density profiles on each stripe can slide relative to 
each other without an energy cost. In other words, there is no rigidity to 
shear deformations of the charge configuration on nearby stripes. This is the 
smectic metal phase [115] . a sliding Luttinger liquid. pOQ] 

The fixed point action for a generic smectic metal phase thus has the form 
(in Fourier space) 

= E^{^+^w)^^}i^(Q)i^ (21) 

where Q = {uj, k, k±), and 9 is the field dual to 0. Here k is the momentum 
along the stripe and k± perpendicular to the stripes. The kernels K{Q) and 
v{Q) are analytic functions of Q whose form depends on microscopic details, 
e. g. at weak coupling they are functions of the inter-stripe Fourier trans- 
forms of the forward and backward scattering amplitudes g2{k±) and g4{k±), 
respectively. In practice, up to irrelevant operators, it is sufficient to keep the 
dependence of the kernels only on the transverse momentum k±. Thus, the 
smectic fixed point is characterized by the effective Luttinger parameter and 
velocity (functions), K{k±) and v{k±). Much like the ordinary ID Luttinger 
liquid, this "fixed point" is characterized by power-law decay of correlations 
functions. This effective field theory also yields the correct low energy de- 
scription of the quantum Hall stripe phase of the 2DEG in large magnetic 
fields. [Ml EH EH [ni fT27] 

In the presence of a spin gap, single electron tunneling is irrelevant |57|. and 
the only potentially relevant interactions involving pairs of stripes a, a' are 
singlet pair (Josephson) tunneling, and the coupling between the CDW order 
parameters. These interactions have the form Hint — J2n (^sc + ^cdw) 
a' — a = n, where 

Sa+n)] 

<l>a+n)]. (22) 

Here Jn are the inter-stripe Josephson couplings (SC), V„ are the 2kp compo- 
nent of the inter-stripe density-density (CDW) interactions, and A is an ultra- 
violet cutoff, yl ~ 1/a where a is a lattice constant. A straightforward calcu- 
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Fig. 16 Phase diagram for a stripe state with a spin gap. 



lation, yields the scaling dimensions Ai^n = ^sc,n and Z\_i.„ = Z\cDW.ri of 

^sc ^'^'^ ^CDW- 

/\±i,„ = f ^ [<ki_)f' (1 - cosnfci) , (23) 

where K{k±) = K{0,0,k±). Since K(fcj_) is a periodic function of k± with 
period 27r, K(fcj_) has a convergent Fourier expansion of the form K{k±) = 
'^iT' cos nk±. We will parametrize the fixed point theory by the coefficients 
Km which are smooth non- universal functions. In what follows we shall discuss 
the behavior of the simplified model with K{k±) = + KiCosfcj_. Here, kq 
can be thought of as the intra-stripe inverse Luttinger parameter, and ni is a 
measure of the nearest neighbor inter-stripe coupling. For stability we require 

Since it is unphysical to consider longer range interactions in Hint than 
are present in the fixed point Hamiltonian, we treat only perturbations with 
n = 1, whose dimensions arc 



Ascs = Asc = Ko-Y (24) 



and 



^CDW4 = -^CDW = -7 N ■ (25) 

For a more general function K{k±), operators with larger n must also be 
considered, but the results are qualitatively unchanged [1641 1209j .^^ 



Asc,2 is the most relevant operator. For a model with K{k±) = [kq + ki cos{k ±)]'^ , 
all perturbations are irrelevant for large kq and small |ko ~ ^^il- 
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In Fig. [in] we present the phase diagram of this model. The dark AB curve 
is the set of points where Aqdw = ^Sd a-nd it is a hne of first order tran- 
sitions. To the right of this hne the inter-stripe CDW couphng is the most 
relevant perturbation, indicating an instability of the system to the forma- 
tion of a 2D stripe crvstal [Tl5] . To the left, Josephson tunneling (which still 
preserves the smectic symmetry) is the most relevant, so this phase is a 2D 
smectic superconductor. (Here we have neglected the possibility of coexis- 
tence since a first order transition seems more likely). Note that there is a 
region of kq ^ Ij and large enough ki, where the global order is supercon- 
ducting although, in the absence of inter-stripe interactions (which roughly 
corresponds to ki =0), the superconducting fluctuations are subdominant. 
There is also a (strong coupling) regime above the curve CB where both 
Josephson tunneling and the CDW coupling are irrelevant at low energies. 
Thus, in this regime the smectic metal state is stable. This phase is a 2D 
smectic non- Fermi liquid in which there is coherent transport only along the 
stripes. 

To go beyond this description we need to construct an effective theory of 
the two-dimensional ordered phase. For instance, the superconducting state 
is a 2D striped superconductor, whereas the crystal is a bidirectional charge 
density wave. A theory of these 2D ordered phases can be developed by 
combining the quasi-one-dimensional renormalization group with an effective 
inter-stripe mean field theory, as in Ref. |26| . which in turn can be fed into a 
2D renormalization group theory. [2] One advantage of this approach is that 
the inter-stripe mean field theory has the same analytic structure as the 
dimensional crossover RG (see Ref. [5]). 

Let us consider the superconducting state, a striped superconductor. In the 
way we constructed this state all ladders are equivalent. Hence this is a period 
2 stripe (columnar) SC phase, similar to the one discussed by Voita. pTO] 
Let us use inter-stripe mean field theory to estimate the critical tempera- 
ture of the 2D state. For the isolated ladder, Tc = as required by the 
Mermin- Wagner theorem. If the inter-stripe Josephson and CDW couplings 
are non-zero, ^ and V 7^ 0, the system will now have a finite SC critical 
temperature, Tc > 0. Now, for x < 0.1, CDW couplings are irrelevant as in 
this range 1/2 < /v < 1. Hence, in the same range, the inter- ladder Josephson 
coupling are relevant and lead to a SC state in a small x with a somewhat 
low Tc which, in inter-stripe (or 'chain') mean field theory can be estimated 

by 

2Jxsc{Tc) = 1 (26) 

In this regime, however, Tc oc St x and it is low due to the low carrier density. 
Conversely, for larger x, K > 1 and xcdw is more strongly divergent than 
Xsc- Thus, for X > 0.1 the CDW couplings become more relevant. This leads 
to an insulating incommensurate CDW state with ordering wave number 

P = 2TTX. 



Electronic Liquid Crystal Phases in Strongly Correlated Systems 



29 









t 


t 

e 


e 


St 


t' 


St 















- -O- 

t 







t' 

■0 



St 







B 



Fig. 17 Model of a period 4 stripe phase 



In the scenario wc just outlined[6j I114| in the 2D regime the system has a 
first order transition from a superconducting state to a non-superconducting 
phase with charge order. However at large enough inter-stripe forward scat- 
tering interactions both couphngs become irrelevant and there is a quantum 
bicritical point separating both phases from a smectic metal (as depicted 
in FigfTG]). However, an alternative possibility is that instead of a bicritical 
point, we may have a quantum tetracritical point and a phase in which SC 
and CDW orders coexist. 



4 Is Inhomogeneity Good or Bad for Superconductivity? 

The analysis we just did raises the question of whether stripe order (that is, 
some form of spatial charge inhomogeneity) is good or bad for superconduc- 
tivity. This question was addressed in some detail in Refs.jSl 1114] where it 
was concluded that a) there is an optimal degree of inhomogeneity at which 
Tc reaches a maximum and b) that charge order in a system with a spin 
gap can provide a mechanism of "high temperature superconductivity" (the 
meaning of which we will specify below) . 

The argument goes as follows. Consider a system with a period 4 stripe 
phase, consisting of an alternating array of inequivalent A and B type ladders 
in the Luther-Emery regime. "'^^ The inter-stripe mean field theory estimate for 



In Ref. |71| a similar pattern was also considered except that the (say) stripes 
do no have a spin gap. This patterns was used to show how a crude model with nodal 
quasiparticles can arise in an inhomogeneous state. 
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the superconducting and CDW critical temperatures now takes the somewhat 
more complex form: 

i2J)\iciT,)xiciTc) = l (27) 
for the superconducting Tc, and 

m\^^^{P,T,)xEDwiP,Tc) = l (28) 

for the CDW Tc- In particular, the 2D CDW order is greatly suppressed 
due to the mismatch between ordering vectors, Pa and Pg, on neighboring 
ladders 

For inequivalent A and B ladders SC beats CDW if the corresponding 
Luttinger parameters satisfy the inequalities 



Ka 



2 > K] 



Kf 



(29) 



The SC critical temperature is then found to obey a power law scaling form 
(instead of the essential singularity of the BCS theory of superconductivity): 



A,. 



2KaKb 



[4KaKb -Ka- Ki 



(30) 



A simple estimate of the effective inter-stripe Josephson coupling, J ^ 5t^ / J 
and of the high energy scale W ^ J , implies that the superconducting critical 
temperature Tc is (power law) small for small with an exponent that 
typically is a ^ 1. 
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Fig. 18 Evolution of the superconducting critical temperature with doping 



These arguments can be used to sketch a phase diagram of the type pre- 
sented in FigU] which shows the qualitative dependence of the SC Tc with 
doping X. The broken line shown is the spin gap As{x) as a function of doping 
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X and, within this analysis, it must be an upper bound on Tj,. Our arguments 
then showed that a period 4 structure can have a substantially larger than 
a period 2 stripe. Consequently, the critical dopings, Xc{2) and a;c(4), for the 
SC-CDW quantum phase transition must move to higher values of x for pe- 
riod 4 compared with period 2. On the other hand, for x ^ Xc the isolated 
ladders do not have a spin gap, and this strong coupling mechanism is no 
longer operative. 

How reliable are these estimates? What we have are mean-field estimates 
for Tc and it is an upper bound to the actual Tc. As it is usually the case, 
Tc should be suppressed by phase fluctuations by up to a factor of 2. On the 
other hand, perturbative RG studies for small J yield the^ame power law 
dependence. This result is asymptotically exact for J « W. Since is a 
smooth function of St/ J', it is reasonable to extrapolate for 5t ~ JT. Hence, 
jrmax ^ g^j-^jj have a "high Tc" . This is in contrast to the exponentially 
small Tc obtained in a BCS-likc mechanism. 

Now, having convinced ourselves that a period 4 stripe will have a larger 
SC Tc than a period 2 stripe one may wonder if an even longer period stripe 
state would do better. It is easy to see that there will be a problem with this 
proposal. Clearly, although the argument we just presented would suggest 
that the exponents will also be of order 1 for longer periods, the problem 
now is that the effective couplings become very small very quickly as the 
Josephson coupling has an exponential dependence on distance (tunneling!) . 
Thus, there must be an optimal period for this mechanism and it is likely to 
be a number larger than 2 but smaller than (say) 6. 

In summary, we have shown that in systems with strong repulsive in- 
teractions (and without attractive interactions), an (inhomogeneous) stripe- 
ordered state can support a 2D superconducting state with a high critical 
temperature, in the sense that it is not exponentially suppressed, with a high 
paring scale (the spin gap). This state is an inhomogeneous version of the 
RVB mechanism. [H 11131 1130] The arguments suggest that there is an opti- 
mal degree of inhomogeneity. There is suggestive evidence in ARPES data in 
La2-a;Baa;Cu04 that show a large pairing scale in the stripe- ordered state 
which support this picture. [204| [85] 



5 The Striped Superconductor: a Pair Density Wave 
state 

We now turn to a novel type of striped superconductor, the pair density wave 
state. Berg et al^5\ \17\ [18] have recently proposed this state as a symmetry- 
based explanation of the spectacular dynamical layer decoupling seen in 
stripe-ordered La2-2;Ba:rCu04 (and Lai rNdn ztSr,rCuOzt) |1331 11991 
and in La2-a;Sr2;Cu04 in magnetic fields. |187| 

Summary of experimental facts for La2-a;Ba2;Cu04 near x = 1/8; 
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• ARPES finds an anti-nodal d-wave SC gap that is large and unsupprcssed 
at a; = 1/8. Hence, there is a large pairing scale in the stripe-ordered state. 

• Resonant X-Ray scattering finds static charge stripe order for T < 

^charge — 54_/^ . 

• Neutron Scattering finds static Stripe Spin order T < Tspin ~ A2K. 

• The in-plane resistivity pab drops rapidly to zero from Tspin to Tkt (the 
Kosterlitz-Thouless (KT) transition). 

• Pab shows KT behavior for Tspin > T > Tkt- 

• pc increases as T decreases for T > T** « 35K. 

• Pc — >■ as T — )- T^D = lOK (the bulk 3D resistive transition). 

• Pc/Pab oo for Tkt >T > T^d- 

• Theres is a Meissner state only below = AK. 

How do we understand these remarkable effects that can be summarized as 
follows: There is a broad temperature range, T^d <T < T2D with 2D super- 
conductivity but not in 3D, as if there is not interlayer Josephson coupling. 
In this regime there is both striped charge and spin order. This can only 
happen if there is a special symmetry of the superconductor in the striped 
state that leads to an almost complete cancellation of the c-axis Josephson 
coupling. 

What else do we know? The stripe state in the LTT ("low temperature 
tetragonal") crystal structure of La2-2;Bai,Cu04 has two planes in the unit 
cell. Stripes in the 2nd neighbor planes are shifted by half a period to minimize 
the Coulomb interaction: 4 planes per unit cell. The anti-ferromagnetic spin 
order suffers a tt phase shift accross the charge stripe which has period 4. Berg 
et aZ[Tn] proposed that the superconducting order is also striped and also 
suffers a tt phase shift. The superconductivity resides in the spin gap regions 
and there is a tt phase shift in the SC order across the anti-fcrromagnctic 
regions. 

The PDW SC state has intertwined striped charge, spin and supercon- 
ducting orders. 

How does this state solve the puzzle? If this order is perfect, the Josephson 
coupling between neighboring planes cancels exactly due to the symmetry 
of the periodic array of tt textures, i.e. the spatial average of the SC order 
parameter is exactly zero. The Josephson couplings Ji and J2 between planes 
two and three layers apart also cancel by symmetry. The first non- vanishing 
coupling J3 occurs at four spacings. It is quite small and it is responsible for 
the non-zero but very low T,. Defects and/or discommensurations gives rise 
to small Josephson coupling Jq neighboring planes. 

Are there other interactions? It is possible to have an inter-plane bi- 
quadratic coupling involving the product SC of the order parameters between 

^® While there is some numerical evidence for a state of this type in variational 
Monte Carlo calculations |87| and in slave particle mean field theorv [l76l I23| (see, 
however, Ref. | 226) and | 138p . a consistent and controlled microscopic theory is yet to 
be developed. Since the difference between the energies of the competing states seen 
numerically is quite small one must conclude that they are all reasonably likely. 
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Fig. 19 Period 4 Striped Superconducting State 



neighboring planes A1A2 and the product of spin stripe order parameters also 
on neighboring planes Mi • M2. However in the LTT structure Mi • M2 = 
and there is no such coupling. In a large enough perpendicular magnetic field 
it is possible (spin flop transition) to induce such a term and hence an effective 
Josephson coupling. Thus in this state there should be a strong suppression 
of the 3D SC Tc but not of the 2D SC T^. 

On the other hand, away from x = 1/S there is no perfect commensura- 
tion. Discommcnsurations arc defects that induce a finite Josephson coupling 
between neighboring planes Ji \x — l/8p, leading to an increase of the 3D 
SC Tc away from x = 1/8. Similar effects arise from disorder which also lead 
to a rise in the 3D SC Tc. 



5.1 Landau- Ginzburg Theory of the Pair Density Wave 

In what follows we will rely heavily on the results of Refs.[l7l [TBI HH]- We 
begin with a description of the order parameters: 

1. PDW (Striped) SC: 

Z\(r) ^ Z\Q(r)e''Q- + Z\_Q(r)e-'Q-'- (31) 

complex charge 2e singlet pair condensate with wave vector Q, (i.e. an 
FFLO type state at zero magnetic field) 

2. Nematic: detects breaking of rotational symmetry: N, a real neutral 
pseudo-scalar order parameter 

3. Charge stripe: pK, unidirectional charge stripe with wave vector K 



A state that is usually described as a pair crystal is commonly known as a pair 
density wave. [H21ll48| However that state cannot be distinguished by symmetry from 
a (two) CDWs coexisting with a uniform SC. 
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4. Spin stripe order parameter: Sq, a neutral complex spin vector order pa- 
rameter. 

These order parameters have the following transformation properties under 
rotations by 7r/2, 

1. The nematic order parameter changes sign: N — > ~N 

2. The CDW ordering wave vector rotates: pK 

3. The SDW ordering wave vector also rotates: Sq Stj^^^q 

4. The striped SC (s or d wave) order parameter: A±q — > ±Z\±7j^^2q (+ for 
s-wave, — for d- wave) 

and by translations by R 

N^N, PK ^ e'^^ '^PK, SQ->e"^-^SQ (32) 

The Landau-Ginzburg free energy functional is, as usual, a sum of terms 
of the form 

J- = J-2 + -7^3 + -^4 + • ■ • (33) 

where J^2i the quadratic term, is simply a sum of decoupled terms for each 
order parameter. There exist a number of trilinear terms mixing several of 
the order parameters described above. They arc 

= Is [p-kSq • Sq + P_kSq • Sq + C.C.] (34) 

+74[p_K/i* Q^iQ + p_K^* Q^Q + C.C.] 

+g^N[A^A^ + Z\* qZ\-q - Z\^Z\q - /i* q/\-q] 

+5.A^[S_Q • Sq - S_Q • Sq] 

where Q = 72,^/2Q, and K ~ 7?.^/2K. The fourth order term, which is more 
or less standard, is shown explicitly below. 

Several consequences follow directly from the form of the trilinear terms, 
Ea. (j35p . One is that, at least in a fully translationally invariant system, the 
first two terms of Eq.pSj) implies a relation between the ordering wave vec- 
tors: K = 2Q. Also, as we will see below, these terms implies the existence 
of vortices of the SC order with half the flux quantum. 

Another important feature of the PDW SC is that it implies the existence 
of a non-zero charge 4e uniform SC state. Indeed, if we denote by the 
(uniform) charge 4e SC order parameter, then the following term in the LG 
expansion is allowed 

(35) 



•^3 = ff4 



Al AoA 



rotation 



by- 

^ 2 



c.c. 
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Hence, the existence of striped SC order (PDW) implies the existence uniform 
charge 4e SC order!^° 

We should also consider a different phase in which there are coexisting 
uniform and striped SC orders, as it presumably happens at low temperatures 
in La2_a;Baa;Cu04. If this is so, there is a non-zero PDW SC order parameter 
Aq as well as an uniform (d-wave) SC order parameter Aq which are coupled 
by new (also trilinear) terms in the LG free energy of the form 

J^3,u = jaA'^{pqA-q + p-qAq) + gpP-2QPQ + ^ rotation + c.c. (36) 

If Z\o 7^ and Aq ^ 0, there is a new pQ component of the charge order!. 
Also, the small uniform component Aq removes the sensitivity to quenched 
disorder of the PDW SC state. 



5.2 Charge 4e SC order and the topological excitations 
of the PDW SC state 

If there is a uniform charge 4e SC order, its vortices must quantized in units 
of hc/Ae instead of the conventional hc/2e flux quantum. Hence, half- vortices 
are natural in this state. To see how they arise let us consider a system deep 
in the PDW SC state so that the magnitude of all the order parameters is 
essentially constant, but their phase may vary. Thus we can write the PDW 
SC order parameter as 

Z\(r) = |Z\q| e'Q-r+*»±QM + \A_^\ e-''Q-r+»e-Q(r) (37) 

where (by inversion symmetry) \Aq\ = \ A^q\ =const. It will be convenient 
to define the new phase fields 0±{r) by 

e±Q{r) = lje+{r)±e^{r)) (38) 

Likewise, in the same regime the CDW order parameter can be written as 

p(r) = IpkI cos(K-r + 0(r)) (39) 

(and a similar expression for the SDW order parameter.) In this notation, 
the second trilinear term shown in Ea. (|35|) takes the form 

-^3.7 = 2^a\p-kAqA^q\ cos(20_(r) - 0(r)) (40) 



A charge 4e SC order parameter is an expectation value of a four fermion operator. 
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Hence, the relative phase 9_ is locked to 0, the Goldstonc boson of the CDW 
(the phason), and are not independently fluctuating fields. Furthermore, the 
phase fields 9±q are defined modulo 2t: while 6^ is defined only modulo tt. 

This analysis implies that the allowed topological excitations of the PDW 
SC are 

1. A conventional SC vortex with Z\6'+ = 2tt and A(j> = 0, with topological 
charges (1,0). 

2. A bound state of a 1/2 vortex and a CDW dislocation, A9^ — tt and 
A(j) — 2tt, with topological charges (±1/2, ±1/2) (any such combination is 
allowed). 

3. A double dislocation, A0+ — and Acj) = An, with topological charge 
(0,1). 

All three topological defects have logarithmic interactions. 

There are now three different pathways to melt the PDW SC[16], de- 
pending which one of these topological excitations becomes relevant (in the 
Kosterlitz-Thouless RG sense |119] ) first. At temperature T, the scaling di- 
mension of a topological excitation of topological charge {p, q) is 

\,q = 7^ {PscP^ + KCDWq^) (41) 

where p^c is the superfiuid density (the stiffness of the 9^ phase field) and 
KcDW is the CDW stiffness (that is, of the phase field). As usual the 
criterion of relevance is that an operator that creates and excitation is relevant 
if its scaling dimension is equal to the space dimension (for details see, for 
instance, Ref.[24]) which in this case is 2. This condition, Ap^q = 2 for each 
one of the topological excitations listed above, leads to the phase thermal 
phase diagram shown in Fig|20l^^ 

Thus, the PDW state may thermally melt in three possible ways: 

1. First into a CDW phase (by proliferating conventional SC vortices, a (1,0) 
topological excitation, followed by a subsequent melting of the CDW into 
the normal (Ising nematic) high temperature phase. This scenario corre- 
sponds to the right side of the phase diagram and , presumably, is what 
happens in La2-a;Baa;Cu04. 

2. A direct melting into the normal (Ising nematic) phase by proliferation of 
fractional vortices, with topological charge (±1/2, ±1/2). 

3. Melting into a charge 4e uniform SC phase by proliferation of double 
dislocations, with topological charge (0, 1). 

The prediction that the PDW state should effectively have a uniform 
charge 4e SC order with an anomalous hc/Ae flux quantization leads to a 
direct test of this state, this can be done by searching for fractional vortices, 

•^^ A more elaborate version of this phase diagram, based on a one-loop Kosterlitz 
RG calculation for for physically very different systems with the same RG structure 
system, was given in Refs. p^TI 11671 ri78 | . 
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Fig. 20 Schematic phase diagram of the thermal melting of the PDW state. 



and similarly of fractional periodicity in the Josephson effect (and Shapiro 
steps). Similarly, the prediction that in the phase in which an uniform (d- 
wave) SC is present there should be a charge-ordered state with period equal 
to that of the SC (and of the SDW) is another direct test of this theory. 



6 Nematic Phases in Fermi Systems 

We now turn to the theory of the nematic phases. The nematic phase is the 
simplest of the liquid crystal states. In this state the system is electronically 
uniform but anisotropic. There are two ways to access this phase. One is by 
a direct transition from the isotropic electronic fluid. The other is by melting 
(thermal or quantum mechanical) the stripe phase. We will consider both 
cases. We will begin with the first scenario in its simplest description as a 
Pomeranchuk instability of a Fermi liquid. 



6.1 The Pomeranchuk Instability 

The central concept of the Landau theory of the Fermi liquid [13] is the quasi- 
particle. A Landau quasiparticle is the elementary excitation of a Fermi liquid 
with the same quantum numbers as a non-interacting electron. A necessary 
condition for the Landau theory to work is the condition that the quasiparti- 
cle becomes sharp (or well defined) at asymptotically low energies, i.e. as the 
Fermi surface is approached. For the Landau quasiparticle to be well defined 



38 



Eduardo Fradkin 



it is necessary that the quasiparticle width, i.e. the quasiparticlc scattering 
rate, to be small on the scale of the quasiparticle energy. The quasiparticle 
scattering rate, the imaginary part of the electron self energy, i7"(w,p), is 
determined by the quasiparticle interactions, which in the Landau theory of 
the Fermi liquid are parametrized by the Landau parameters. Except for the 
BCS channel, the forward scattering interactions (with or without spin flip) 
are the only residual interactions among the quasiparticlcs that survive at 
low energies. |169| I190j 

The Landau "parameters" are actually functions F'^'"^(p, p') quantify the 
strength of the forward scattering interactions among quasiparticlcs at low 
energies with momenta p and p' close to the Fermi surface in the singlet 
(charge) channel (S) or the triplet (spin) channel (A). For a translation- 
ally invariant system it depends only on the difference of the two momenta, 
F^PtV') = P{v ~ P')- Furthermore, if the system is also rotationally invari- 
ant, the Landau parameters can be expressed in an angular momentum basis. 
In 3D they take the form fI'^ (with £ = 0,1,2,... and |m| < i) , while in 2D 
they are simply i^,fj'"^ (where m ^ 1). We will see below that in some cases 
of interest we will also need to keep the dependence on a small momentum 
transfer in the Landau parameters {i.e. p and p' will not be precisely at the 
FS) even though it amounts to keeping a technically irrelevant interaction. 
On the other hand, for a lattice model rotational invariance is always broken 
down to the point (or space) group symmetry of the lattice. In that case the 
Landau parameters are classified according to the irreducible representations 
of the point (or space) group of the lattice, e.g. the C4 group of the square 
lattice. 

It is well known in the Landau theory of the Fermi liquid that the ther- 
modynamic stability of the Fermi liquid state requires that the Landau pa- 
rameters cannot be too negative. This argument, due to Pomeranchuk [T7T] . 
implies that if in one channel the forward scattering interaction becomes suf- 
ficiently negative (attractive) to overcome the stabilizing effects of the Pauli 
pressure, the Fermi liquid becomes unstable to a distortion of the FS with 
the symmetry of the unstable channel. 

Oganesyan et aZ }163j showed that in a 2D system of interacting fermions, 
the Pomeranchuk instability in fact marks a quantum phase transition to 
a nematic Fermi fluid. We will discuss this theory below in some detail. 
While the the theory of Oganesyan and coworkers applied to a system in the 
continuum, Kee and coworkers jl06l ITOTl considered a lattice model. Hints of 
nematic order in specific models had in fact been discovered independently 



Although the Pomeranchuk argument is standard and reproduced in all the text- 
books on Fermi liquid theory (see, e.g. Ref. |13|) the consequences of this instability 
were not pursued until quite recently. 
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(but not recognized as such originally), notably by the work of Metzner and 
coworkers [Ha Ei [JSOl ttSl 12151 HSJ^^ 

There is by now a growing literature on the nematic instability. Typically 
the models, both in the continuum [TBS] or on different lattices [l06[ 11231 11231 
I175j , are solved within a Hartree-Fock type approximation (with all the limi- 
tations that such an approach has), or in special situations such as vicinity to 
Van Hove singularities j 1 7) I225j and certain degenerate band crossings [196] 
(where the theory is better controlled) , or using uncontrolled approximations 
to strong coupling systems such as slave fermion/boson methods. |2231 1154] . 
A strong coupling limit of the Emery model of the cuprates was shown to 
have a nematic state in Ref. |116j (we will review this work below). Finally 
some non-perturbative work on the nematic quantum phase transition has 
been done using higher dimensional bosonization in Refs. fl2(]l 1128] and by 
RG methods. [HI] 

Extensions of these ideas have been applied to the problem of the ne- 
matic phase seen in the metamagnetic bilayer ruthenate Sr3Ru207 rely- 
ing either on the van Hove mechanism [TUS] 12221 1173| I174j or on an orbital 
ordering mechanism jl791 1213j . and in the new iron-based superconducting 
compounds. [60] I221j More recently nematic phases of different types have 
been argued to occur in dipolar Fermi gases of ultra-cold atoms. [571 166p "^ 



6.2 The Nematic Fermi Fluid 



Here I will follow the work of Oganesyan, Kivelson and Fradkin [163j and 
consider first the instability in the charge (symmetric) channel. Oganesyan 
et al defined a charge nematic order parameter for two-dimensional Fermi 
fluid is the 2x2 symmetric traceless tensor of the form 

Q(.).-j^*tw(f-^.'|:^..)*H (42) 

It can also be represented by a complex valued field Q2 (x) whose expectation 
value is the nematic phase is 

(Q2) = ^ (a, + idyf ^) = IQ2I e"''^ = Qn + tQi2 ^ (43) 

Q2 transforms under rotations in the representation of angular momentum 
£ = 2. Oganesyan et al showed that if (Q2) 7^ then the Fermi surface 

In fact, perturbative renormalization group calculations |93 1 194 | have found a run- 
away flow in the d^2_y2 particle-hole channel, which is a nematic instability, but it 
was not recognize it as such. See, however, Ref.|84|. 

Another class of nematic state can occur inside a d^2_y2 superconductor. This 
quantum phase transition involves primarily the nodal quasiparticles of the super- 
conductor and it is tractable within large type approximations. | 108l I99| 
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spontaneously distorts and becomes an ellipse with eccentricity cx Q. This 
state breaks rotational invariance mod tt. 

More complex forms of order can be considered by looking at particle-hole 
condensates with angular momenta £ > 2 (see Ref. |194j ) 

{Qe) = {^Hd. + ^dyY^) (44) 

For £ odd, this condensate breaks rotational invariance (mod 2tt/£). It also 
breaks parity V and time reversal T but VT is invariant. For example 
the condensate with ^ = 3 is effectively equivalent to the "Varma loop 
state" . [2051 1206] . The states with £ even are also interesting, e.g. a hexatic 
state is described by a particle- hole condensate with £ = 6.[T0] 

In a 3D system, the anisotropic state is described by an order parameter 
Qij which is a traceless symmetric tensor (as in conventional liquid crystals 
H5] ) . More generally, we can define an order parameter that transforms 
under the {£, m) representation of the group of 5*0(3) spatial rotations. 

Oganesyan et al considered in detail Fermi liquid type model of the nematic 
transition and developed a (Landau) theory of the transition ("Landau on 
Landau" ) . The Hamiltonian of this model describes (spinless) fermions in the 
continuum with a two-body interaction corresponding to the £ — 2 particle- 
hole angular momentum channel. The Hamiltonian is 

H = j dv If ^(r)e(v)*(r) + J / dr j dv'¥^{v - r')Tr[Q(r)Q(r')] (45) 

where the free-fermion dispersion (near the FS) is e(k) = vpql^ + o-i'^^)'^] 
(here q = |k| — kp), and the interaction is given in terms of the coupling 

F^{v) = (27r)-2 / dcie^^-^ \ ^ (46) 

where F2 is the £ ~ 2 Landau parameter, and k measures the range of these 
interactions. Notice that we have kept a cubic momentum dependence in the 
dispersion, which is strongly irrelevant in the Landau Fermi liquid phase (but 
it is needed to insure stability in the broken symmetry state). 

The Landau energy density functional for this model has the form (which 
can be derived by Hartree-Fock methods or, equivalently, using a Hubbard- 
Stratonovich decoupling) 

W[Q] = f (Q) - ^Tr[QDQ] - ^TrlQ^DQ] + . . . (47) 

where k, and k are the two effective Franck constants (see Ref.j29]). The 
uniform part of the energy functional, £(Q), is given by 

8{Q) = f (0) + ^Tr[Q2] + |Tr[Q4] + . . . (48) 
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^=^ + F^ (49) 

Np is the density of states at the Fermi surface, and the coefficient of the 
quartic term is 

3aNF\F2\' 

Ep = vpkp is the Fermi energy. [f63j The (normal) Landau Fermi hquid 
phase is stable provided A > 0, or, what is the same, if 2NpF2 > —I which 
is the Pomeranchuk condition (in this notation) . On the other hand, thermo- 
dynamic stability also requires that B > 0, which implies that the coefficient 
of the cubic correction in the dispersion be positive, a > 0. If this condition 
is not satisfied, as it is the case in simple lattice models [TUS] . higher order 
terms must be kept to insure stability. However, in this case the transition is 
typically first order. 

This model has two phases: 

• an isotropic Fermi liquid phase, ^ > 

• a nematic (non-Fermi liquid) phase, ^ < 

separated by a quantum critical point at the Pomeranchuk value, 2NpF2 — 
-1. 

Let us discuss the quantum critical behavior. We will parametrize the 
distance to the Pomeranchuk QCP by 

and define s ~ u/qvp. The transverse collective nematic modes have Landau 
damping at the QCP. [T55] Their effective action has the form 



S^= f dujdq (nq^ + S ^ i-^) |Qi(w,q)|^ 
J qvp 



(52) 



which implies that the dynamic critical exponent is z = 3.^^ 

According to the standard perturbative criterion of Hertz [55] and Millis jl51j . 
the quantum critical behavior is that of an equivalent type field theory in 
dimensions D ^ d + z which in this case is I? = 5. Since the upper critical 
(total) dimension is 4, the Hertz-Millis analysis would predict that mean field 
theory is asymptotically exact and that the quartic (and higher) powers of 
the order parameter field are irrelevant at the quantum critical point (for an 
extensive discussion see Ref. [186j .) However we will see below that while this 



■^^ There are other collective modes at higher energies. In particular there is an 
underdamped longitudinal collective mode with z = 2. |163| These higher energy 
modes contribute to various crossover effects [25IJ] . but decouple in the asymptotic 
quantum critical regime. 
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analysis is correct for the bosonic sector of the theory, i.e. the behavior of 
the bosonic collective modes such as the order parameter itself, the situation 
is far less clear in the fermionic sector. We will come back to this question 
below. 

Let us discuss now the physics of the nematic phase. In the nematic phase 
the FS is spontaneously distorted along the direction of the (director) order 
parameter (see Fig l21|) and exhibits a quadrupolar (d-wave) pattern, i.e. the 
Fermi wave vector has an angular dependence kpiO) oc cos 29 (in 2D). Indeed, 
in the nematic phase the Hartree-Fock wave function is a Slater determinant 
whose variational parameters determine the shape of the FS. 




Fig. 21 Spontaneous distortion of the Fermi surface in the nematic phase of a 2D 
Fermi fluid. 

In principle, a wave function with a similar structure can be used to suggest 
(as it was done in Ref. |163| ) that it should also apply to the theory of the 
electronic nematic state observed in the 2DEG in large magnetic fields. In 
that framework one thinks of the 2DEG in a half-filled landau level as an 
equivalent system of "composite fermions" |101| , fermions coupled to a Chern- 
Simons gauge field [139[|8T] . It has been argued [l83] that this state can be well 
described by a Slater determinant wave function, projected onto the Landau 
level. The same procedure can be applied to the nematic wave function, and 
some work has been done along this direction [49]. A problem that needs to be 
solved first is the determination of the Landau parameters of the composite 
fermions of which very little (that makes sense) is known. 

A simple (Drude) calculation then shows that the transport is anisotropic. 
The resistivity in the nematic phase, due to scattering from structureless 
isotropic impurities, yields the result that the resistivity tensor is anisotropic 
with an anisotropy controlled by the strength of the nematic order parameter: 
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Pxx Pyy _ 1 my rrix _ Rc Q _^ 0{Q^) (53) 

Pxx + Pyy 2my +mx Ep 

where rrix and m,y are the (anisotropic) effective masses of the quasiparticles 
in the nematic state. In general it is a more complex odd function of the 
order parameter. 

In the nematic phase the transverse Goldstone boson is generically over- 
damped (Landau damping) except for a finite set of symmetry directions, 
(j) = 0,±7r/4,±7r/2, where it is underdamped. Thus, 2 = 3 scaling also applies 
to the nematic phase for general directions. Naturally, in a lattice system the 
rotational symmetry is not continuous and the transverse Goldstone modes 
are gapped. However, the continuum prediction still applies if the lattice sym- 
metry breaking is weak and if either th energy or the temperature is larger 
that the lattice anisotropy scale. 

On the other hand, the behavior of the fermionic correlators is much 
more strongly affected. To one loop order, the quasiparticle scattering rate, 
E (u), p) is found to have the behavior 



4/3 



2vFkF 



2/3 



(54) 



for k along a general direction. On the other hand, along a symmetry direction 



3Nfk {nkiy/"^ 



3/2 



(55) 



Hence, the entire nematic phase is a non-Fermi liquid (again, with the caveat 
on lattice symmetry breaking effects). 

At the Pomeranchuk quantum critical point the quasiparticle scattering 
rate obeys the same (one loop) scaling shown in Eq. (|54| . S (w) cx |cjp/'^, both 
in continuum fl63| and lattice models}^ . but it is isotropic. In the quantum 
critical regime the electrical resistivity obeys a T^/^ law. [47]. Also, both in 
the nematic phase (without lattice anisotropy) and in the quantum critical 
regime, the strong nematic fluctuations yield an electronic contribution to 
the specific heat that scales as T^/^ (consistent with the general scaling form 
T''/^[THS]) which dominates over the standard Fermi liquid linear T depen- 
dence at low temperatures. [13] 

Since S"{uj) ^ ^'(w) (as w 0), we need to asses the validity of these 
results as they signal a failure of perturbation theory. To this end we have used 
higher dimensional bosonization as a non-perturbative tool[79l [27l [28l [Ml 
[95] . Higher dimensional bosonization reproduces the collective modes found 
in Hartree-Fock-|- RPA and is consistent with the Hertz-Millis analysis of 
quantum criticality: des = d + z = 5. |126[ 1128] . Within this approach, the 
fermion propagator takes the form 



GF{x,t) = Go{x,t)Z{x,t) 



(56) 
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At the Nematic-FL QCP it has the form 

Gf{x, 0) = Go(x, 0) e-^onst. \x\'l^ (57) 

at equal times, and 

Gf(0, t) = Go(0, t) e-const. Kp^/^ \nt (gg) 

at equal positions. Notice that these expressions are consistent with the ex- 
pected z = 3 scaling even though the time and space dependence is not a 
power law. The quasiparticle residue is then seen to vanish at the QCP: 

Z = lim Z(.T,0) = (59) 

a;— >oo 

However, the single particle density of states, Niuj) = — ilmG(aj,0), turns 
out to have a milder behavior: 

N{uj) = A^(0) (1 - const'. Inw) (60) 



n{k) 




Let us now turn to the behavior near the QCP. For T = and (5^1 (on 
the Fermi Liquid side) the quasiparticle residue is now finite 

Z oc e-co^st-/^ (61) 

but its dependence on the distance to the nematic QCP is an essential singu- 
larity. On the other hand, right at the QCP ((5 = 0), and for a temperature 
range Tp ^ T ^ T^, the equal-time fermion propagator is found to vanish 
exactly 

Zix, 0) oc e-const. Tx'' In {L/x) ^ g as L ^ 00 (62) 
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but, the equal-position propagator Z{0,t) remains finite in the thermody- 
namic hmit, L — > oo\ This behavior has been dubbed "Local quantum crit- 
icality".^^ On the other hand, irrelevant quartic interactions of strength u 
lead to a renormalization of S that smears the QCP at T > |151] 

S^6(T)^ ^uT In [uT^/^^ (63) 

leading to a milder behavior at equal-times 

0) cx e-'^"""*- where C = 5{T)~^/^ (64) 

These results are far from being universally accepted. Indeed Chubukov 
and coworkers [551 [Ml 1 1 8 1 j have argued that the perturbative non- Fermi liquid 
behavior, E [lo) ^ lu"^/^ ^ which is also found at a ferromagnetic metallic QCP, 
persists to all orders in perturbation theory and can recover the results of 
higher dimensional bosonization only by taking into account the most infrared 
divergent diagrams. The same non-Fermi liquid one-loop perturbative scaling 
has been found in other QCPs such as in the problem of fermions (relativistic 
or not) at finite density coupled to dynamical gauge fields. This problem has 
been discussed in various settings ranging from hot and dense QED and 
QCD[92j [121 [20], to the gaugc-spinon system in RVB approaches to high 
superconductors jl82( [100. 159. 170t[31 ]ll31j to the compressible states of the 
2DEG in large magnetic fields [81]. In all cases these authors have also argued 
that the one-loop scaling persists to all orders. In a recent paper Metlitski 
and Sachdev |149| found a different scaling behavior. 

We end with a brief discussion on the results in lattice models of the 
nematic quantum phase transition. This is important since, with the possible 
exception of the 2DEG in large magnetic fields and in ultra-cold atomic 
systems, all strongly correlated systems of interest have very strong lattice 
effects. The main difference between the results in lattice models and in the 
continuum is that in the former the quantum phase transition (at the mean 
field, Hartree-Fock, level) has a strong tendency to be first order. Although 
fluctuations can soften the quantum transition and turn the system quantum 
critical (as emphasized in Ref. jl02j ). nevertheless there are good reasons for 
the transition to be first order more or less generically. One is that if the 
stabilizing quartic terms are negative {e.g. say due to the band structure), this 
also results, in the case of a lattice system, in a Lifshitz transition at which 
the topology of the FS changes from closed to open. This cannot happen in 
a continuous way. 



A similar behavior was found in the quantum Lifshitz model at its QCP [70] . 
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6.3 Generalizations: Unconventional Magnetism and 
Time Reversal Symmetry Breaking 

We will now consider briefly the extension of these ideas to the spin-triplet 
channel j220[ I219| . In addition to particle- hole condensates in the singlet 
(charge) channel we will be interested particle-hole condensates in the spin 
(triplet) channel. In 2D the order parameters for particle-hole condensates in 
the spin triplet channel are (here a, (3 =t, \) 

Qiir) = {^i{r)cjlp (a. + idyf^p{r)) ^ + (65) 

These order parameters transform under both S0{2) spatial rotations and 
under the internal SU{2) symmetry of spin, li £ ^ the state has a broken ro- 
tational invariance in space and in spin space. These states arc a particle-hole 
condensate analog of the unconventional superconductors and supcrfluids, 
such as HC3A and HcsB. Indeed one may call these states "unconventional 
magnetism" as the £ = (isotropic) state is just a ferromagnet. In 2D these 
states are then given in terms of two order parameters, each in the vector 
(adjoint) representation of the SU{2) spin symmetry.^^ We will discuss only 
the 2D case. The order parameters obey the following transformation laws: 

1. Time Reversal: 

rQ?r-i = i-iY+'Qi (66) 

2. Parity: 

VQ^eV-' = i-lYQI (67) 

3. Q1 rotates under an S'Ospin(3) transformation, and transforms as Q1 — )■ 
Qa^iie under a rotation in space by an angle 9. 

4. Q° is invariant under a rotation by tt/^ followed by a spin flip. 

Wu and collaborators |22d[ 1219] have shown that these phases can also be 
accessed by a Pomeranchuk instability in the spin (triplet) channel.^* They 
showed that the Landau-Ginzburg free energy takes the simple form 

F[n] =r(|ni|2 + |n2n+t;i(|ni|2 + |n2n2+t;2|ni xn2p (68) 

The Pomeranchuk instability occurs at r = 0, i.e. for NpFf- = -2 (with i > 
1), where arc the Landau parameters in the spin- triplet channel. Notice 
that this free energy is invariant only by global 5*0(3) rotations involving 
both vector order parameters, ni and 112. Although at this level the 5*0(3) 
invariance is seemingly an internal symmetry, there are gradient terms that 
lock the internal 50(3) spin rotations to the "orbital" spatial rotations. (See 
Ref. |219j ). A similar situation also occurs in classical liquid crystals. [13] 

In 3D the situation is more complex and the possible are more subtle. In particular, 
in 3D there are three vector order parameters involved. | 219 | 

■^^ The 1 = case is, of course, just the conventional Stoner ferromagnetic instability. 
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I = 1 a-phase I = 2 



Fig. 22 The a-phases in the £ = 1 and £ = 2 spin triplet channels. The Fermi 
surfaces exhibit the p and d-wave distortions, respectively. 

At the level of the Landau-Ginzburg theory the system has two phases 
with broken 5*0(3) invariance: 

1. Iiv2 > 0, then the two 5*0(3) spin vector order parameters must be parallel 
to each other, rii x n2 = 0. They dubbed this the "a" phase. In the a phases 
the up and down Fermi surfaces are distorted (with a pattern determined 
by £) but are rotated from each other by tt/£. One case of special interest 
is the a phase with £ = 2. This is the "ncmatic-spin-nematic" discussed 
briefly in Ref. |117) .^^ In this phase the spin up and spin down FS have 
an £ = 2 quadrupolar (ncmatic) distortion but are rotated by 7r/2 (see 
FiglMl). 

2. Conversely, if V2 < 0, then the two S'0(3) spin vector order parameters 
must be orthogonal to each other, rii • n2 = and |ni| = |n2|. Wu et al 
dubbed these the "/?" phases. In the /3 phases there are two isotropic FS 
but spin is not a good quantum number. In fact, the electron self energy in 
the /3-phases acquires a spin-orbit type form with a strength given by the 
magnitude of the order parameter. The mean-field electronic structure thus 
resembles that of a system with a strong and tunable spin-orbit coupling 
(i.e. not of 0{{vf / c)^) as it is normally the case.) 

We can define now a d vector: 

d(k) = (cos(^6'k), sin(£6lk), 0) (69) 

In the P phases the d vector winds about the undistorted FS. For the 
special case oit = 1, the windings arc w ^ 1 (corresponding to a "Rashba" 
type state) and w = ~1 (corresponding to a "Dresselhaus" type state). 
For the d-wave case, the winding numbers are w = ±2 (see Fig. 1231 ) 



The term "nematic-spin- nematic" is a poor terminology. A spin nematic is a state 
with a magnetic order parameter that is a traceless symmetric tensor, which this state 
does not. 



48 



Eduardo Fradkin 




Fig. 23 The /3-phases in tlie 1 = 2 triplet channel. Spin configurations exhibit the 
vortex structure with winding number w = ±2. These two configurations can be 
transformed to each other by performing a rotation around the x-axis with the angle 

of TT. 



These phases have a rich phenomenology of collective modes and topological 
excitations which we will not elaborate on here. See Ref. |219j for a detailed 
discussion.'^" 

Fermionic systems with dipolar magnetic interactions may be a good can- 
didate for phases similar to the ones we just described (See Refs.[67( l66]. and 
it is quite possible that these systems may be realized in ultra-cold atomic 
gases. In that context the anisotropic form of the dipolar interaction provides 
for a mechanism to access some of this physics. Indeed, in the case of a fully 
polarized (3D) dipolar Fermi gas, the FS will have an uniaxial distortion. If 
the polarization is spontaneous (in the absence of a polarizing external field) 
this phase is actually a ferro-nematic state, a state with coexisting ferromag- 
netism and nematic order. If the system is partially polarized then the phase 
is a mix of nematic order and ferromagnetism coexisting with a phase with 
a non-trivial "spin texture" in momentum space. 

It turns out that generalizations of the Pomeranchuk picture of the ne- 
matic state to multi-band electronic systems can describe metallic states 
with a spontaneous breaking of time reversal invariance. This was done in 
Ref. |194] where it was shown that in a two-band system (i.e. a system with 
two Fermi surfaces) it is possible to have a metallic state which breaks time 
reversal invariance and exhibits a spontaneous anomalous Hall effect. While 
the treatment of this problem has a superficial formal similarity with the 
triplet (spin) case, i.e. regarding the band index as a "pseudo- spin" (or fia- 
vor), the physics differs considerably. At the free fermion level the fcrmion 
number on each band is separately conserved, leading to a formal SU{2) 
symmetry. However, the interacting system has either a smaller U{1) x U{1) 



The p-wave {£ = 1) /3 phase has the same physics as the 'spin-split' metal of 
Ref. |91|. A similar state was proposed in Ref. | 207] as an explanation of the "hidden 
order" phase of URu2Si2. 
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invariance or, more generally, Z2 x Z2 invariance, as the more general interac- 
tions preserve only the parity of the band fermion number. 1 194] At any rate it 
turns out that analogs of the "a" and "/?" phases of the triplet channel exist 
in multi-band systems. The "a" phases break time reversal and parity (but 
not the product). An example of such metallic (gapless) states is the "Varma 
loop state" . [2061 1191] ■ The "/3" states break time reversal invariance (and 
chirality). In the "/?" phases there is a spontaneous anomalous Hall effect, 
i.e. a zero field Hall effect with a Hall conductivity that is not quantized as 
this state is a metal'^^, whereas the "a" phases there is not. 



7 Nematic Order in the Strong Correlation Regime 
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Fig. 24 The Emery model of the CuO lattice. 

We will now discuss how a nematic state arises as the exact ground state 
in the strong coupling limit of the Emery model. |116| The Emery model is a 
simplified microscopic model of the important electronic degrees of freedom 
of the copper oxides. [52j In this model, the CuO plane is described as a square 
lattice with the Cu atoms residing on the sites and the O atoms on the links 
(the medial lattice of the square lattice). On each site of the square lattice 
there is a single d^2_y2 Cu orbital, and a (py) O orbital on each site of 
the medial along the x (y) direction. We will denote by dj.(r) the operator 
the creates a hole on the Cu site r and by pj,.^(r -I- ^) and p.y.cr(r + ^) the 
operators the create a hole on the O sites r + ^f■ and r + ^ respectively. 

The Hamiltonian of the Emery model is the sum of kinetic energy and 
interaction terms. The kinetic energy terms consist of the hopping of a hole 
from a Cu site to its nearest O sites (with amplitude tpd), an on-site energy 
e > on the O sites (accounting for the difference in "affinity" between Cu 
and O), and a (small) direct hopping between nearest-neighboring O sites, 
tpp. The interaction terms are just the on-site Hubbard repulsion Ud (on the 



This is consistent with the general arguments of Ref.|80|. 
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Cu sites) and Up (on the O sites) as well as nearest neighbor ("Coulomb") 
repulsive interactions of strength Vpd (between Cu and O) and Vpp (between 
two nearest O). (See FiglMl) It is commonly believed that this model is 
equivalent to its simpler cousin, the one band Hubbard model. However, 
while this equivalency is approximately correct in the weak coupling limit, it 
is known to fail already at moderate couplings. We will see that in the strong 
coupling limit, no such reduction (to a "Zhang- Rice singlet") is possible. 

Let us look at the energetics of the 2D Cu O model in the strong coupling 
limit. By strong coupling we will mean the regime in which the following 
inequalities hold: 

^,^,|^,^^0, U,>Up»Vpa>Vpp, and ^^0 

Up Ud Vpd Vpp tpd 

(70) 

as a function of hole doping a; > 0, where x is the number of doped holes per 
Cu. In this regime, neither Cu nor O sites can be doubly occupied. At half 
filling, X = 0, the holes occupy all the Cu sites and all O sites are empty. 
At half-filling and in this strong coupling regime the Emery model (much as 
the Hubbard model) is equivalent to a quantum Heisenberg antiferromagnet 

St" 

with a small exchange coupling Jh ~ Trrr^- This is the double-exchange 

mechanism. (It turns out that in this model the four-spin ring exchange in- 
teractions can be of the same order of magnitude as the Heisenberg exchange 

jH.mi-) 

Let us consider now the very low doping regime, a; — ^ 0. Any additional 
hole will have to be on an O site. The energy to add one hole {i.e. the chemical 
potential fi of a hole) is /i = 2Vpd + e. Similarly, the energy of two holes on 
nearby O sites is 2/i + Vpp. It turns out that in this strong coupling regime, 
with tpp = 0, the dynamic of the doped holes is strongly constrained and 
effectively becomes one-dimensional. The simplest allowed move for a hole 
on an O site, which takes two steps, is shown in Figi25l The final and initial 



=> 



Fig. 25 An allowed two-step move. 



states are degenerate, and their energy is Eq + fj., where Eq is the ground 
state energy of the undoped system. If this was the only allowed process, the 
system would behave as a collection of ID fermionic systems. 
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To assess if this is correct let us examine other processes to the same 
(lowest) order in perturbation theory (in powers of the kinetic energy). One 
possibility is a process in which in the final state the hole "turned the cor- 
ner" (went from being on an x oxygen to a near y oxygen) . However for that 
to happen it will have to go through an intermediate state such as the one 
shown in Fig|26k. This intermediate state has an energy Eq + ii + Vpp. Hence, 

the effective hopping matrix element to turn the corner is tes = ^ tpd, 
which is strongly suppressed by the Coulomb effects of Vpp. In contrast, the 
intermediate state for the hole to continue on the same row (see Fig lM]) is 

Eq + n + e. Thus the effective hopping amplitude instead becomes toff = 
which is not suppressed by Coulomb effects of Vpp. All sorts of other processes 
have large energy denominators and are similarly suppressed (for a detailed 
analysis see Ref. |116j .) The upshot of this analysis is that in the strong cou- 



Fig. 26 Intermediate states for processes in which a hole a) turns a corner and b) 
continues on the same row. 

pling limit the oped holes behave like a set of one-dimensional fermions (one 
per row and column). The argument for one-dimensional effective dynamics 
can in fact be made more rigorous. In Ref. |116] it is shown that to leading 
order in the strong coupling expansion the system is is a generalized t ~ J 
type model (with J effectively set to zero). In this limit there are an infinite 
number of conserved charges, exactly one per row and one per column. 

The next step is to inquire if, at fixed but very small doping a; — > 0, the 
rows and columns are equally populated or not. Consider then two cases; a) 
all rows and columns have the same fermion density, and b) the columns (or 
the rows) are empty. Case a) is isotropic while case b) is nematic. It turns out 
that due to the effects of the repulsive Coulomb interaction V^p, the nematic 
configuration has lower energy at low enough doping. The argument is as 
follows. For the nematic state (in which all rows are equally populated but 
all rows are empty), the ground state energy has an expansion in powers of 
the doping x of the form: 



i^nematic ^ E{x = 0) + X + W + 0{x^) 



(71) 
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where Ac ~ 2Vpd + e + . . . and W = tt'^ / &m* . The energy of the isotropic 
state (at the same doping x) is 

i^isotropic = E{x = 0) + Zic a; + (1/4) + Feff + . . . (72) 

where Vcr oc Vpp is an effective coupling for holes on intersecting rows and 
columns. Clearly, for a; small enough -^isotropic > ^'ncmatic- Therefore, at low 
enough doping the ground state of the Emery model in the strong coupling 
limit is a ncmatic, in the sense that it breaks spontaneously the point sym- 
metry group of the square lattice, C4.^^ What happens for larger values of 
X is presently not known. Presumably a complex set of stripe phases exist. 
How this analysis is modified by the spin degrees of freedom is an open im- 
portant problem which may illuminate our understanding of the physics of 
high temperature superconductors. 

The nematic state we have found to be the exact ground state is actually 
maximally nematic: the nematic order parameter is 1. To obtain this result 
we relied on the fact that we have set tpp = as this is by far the smallest 
energy scale. The state that we have found is reminiscent of the nematic state 
found in 2D mean field theories of the Pomcranchuk transition [106| in which 
it is found that the nematic has an open Fermi surface, as we have also found. 
Presumably, for the more physical case of tpp ^ 0, the strong ID-like nematic 
state we found will show a crossover to a 2D (Ising) Nematic Fermi liquid 
state. 



8 The Quantum Nematic-Smectic QCP and the Mehing 
of the Stripe Phase 

We will now turn to the problem of the the quantum phase transition between 
electron stripe and nematic phases. For simplicity we will consider only the 
simpler case of the charge stripe and the charge nematic, and we will not dis- 
cuss here the relation with antiferromagnetic stripes and superconductivity. 
Even this simpler problem is not well understood at present. 

In classical liquid crystals there are two well established ways to describe 
this transition, known as the smcctic A-ncmatic transition. One approach is 
the McMillan-de Gennes theory, and it is a generalization of the Landau- 
Ginzburg theory of phase transitions to this problem (sec Ref.|45|.) The 
other approach regards this phase transition as a melting of the smectic 
by proliferation of its topological excitations, the dislocations of the smectic 
order. [1601 [T981 [29] 



It also turns out that in the (so far physically unrealizable) case oi x = 1, the 
ground state is a nematic insulator as each row is now full. However, for a; — > 1 the 
ground state is again a nematic metal. 
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There are however important (and profound) differences between the prob- 
lem of the quantum melting of a stripe phase into a quantum nematic and its 
classical smectic/ncmatic counterpart. The classical problem refers to three- 
dimensional liquid crystals whereas here we are interested in a two- dimen- 
sional quantum system. One may think that the time coordinate, i.e. the 
existence of a time evolution dictated by quantum mechanics, provides for 
the third dimension and that two problems may indeed be closely related. 
Although to a large extent this is correct (this is apparent in the imaginary 
time form of the path integral representation) some important details are 
different. 

The problem we are interested in involves a metal and hence has dynamical 
fermionic degrees of freedom which do not have a counterpart in its classical 
cousin. One could develop a theory of quantum melting of the stripe phase by 
a defect proliferation mechanism by considering only the collective modes of 
the stripe. Theories of this type can be found in Refs. |2251[55] and in Ref. |178j 
(in the context of a system of cold atoms) which lead to several interesting 
predictions. Theories of this type would describe correctly an insulating stripe 
state in which the fermionic degrees of freedom are gapped and hence not part 
of the low energy physics. The problem of how to develop a non-perturbative 
theory of this transition with dynamical fermions is an open and unsolved 
problem. '^^ 

Another important difference is that in most cases of interest the quan- 
tum version of this transition takes place in a lattice system, thus, even if 
the stripe state may be incommensurate, and hence to a first approximation 
be allowed to "slide" over the lattice background, there is no continuous ro- 
tational invariance but only a point group symmetry leftover. Thus, at least 
at the lowest energies, the nematic Goldstone mode which plays a key role 
in the classical problem, is gapped and drops out of the critical behavior. 
However one should not be a purist about this issue as there may be signifi- 
cant crossovers that become observable at low frequencies and temperatures 
if the lattice effects are weak enough. Thus it is meaningful to consider a 
system in which the lattice symmetry breaking are ignored at the beginning 
and considered afterwards. 

In Ref. |195j a theory of the quantum melting of a charge stripe phase is 
developed using an analogy with the McMillan-deGennes theory. The main 
(and important) difference is the role that the fermionic degrees of freedom 
play in the dynamics. Thus we will describe the stripe (which at this level is 
equivalent to a CDW) by its order parameter, the complex scalar field «?(r, t), 
representing the Fourier component of the charge density operator near the 
ordering wavevector Q.'^'* 

We will assume that the phase transition occurs in a regime in which the 
nematic order is well developed and has a large amplitude \Af\. In this regime 



Important work with a similar approach has been done on the problem of the 
quantum melting of the stripe state in quantum Hall systems. | 214l [T77) . 
For a different perspective see Ref. 1 109] . 
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the fluctuations of the amphtude of the nematic order parameter are part 
of the high energy physics, and can be integrated out. In contrast we wiU 
assume that the phase mode of the nematic order, the Goldstone mode, is 
either gapless (as in a continuum system) or of low enough energy that it 
matters to the physics (as if the lattice symmetry breaking is sufficiently 
weak). In this case we will only need to consider the nematic Goldstone (or 
'pseudo- Goldstone') field which we will denote by ip{r,t). 

We should note that there is another way to think about this problem, 
in which one considers the competition between the CDW order parameters 
(two in this case), the nematic order and the normal Fermi liquid near a 
suitable multi-critical point. This problem was considered briefly in Ref. [195| 
and revisited in more detail in Ref. |153j . The main conclusion is that for a 
(square) lattice system, the multicritical point is inaccessible as it is replaced 
by a direct (and weak) fluctuation-induced first-order transition from the FL 
to the CDW state. Thus, the theory that we discuss here applies far from 
this putative multicritical point in a regime in which, as we stated above, the 
nematic order is already well developed and large. 

Following Ref. |195j we will think of this quantum phase transition in the 
spirit of a Hertz-Millis type approach and postulate the existence of an order 
parameter theory coupled to the fermionic degrees of freedom. The quan- 
tum mechanical action of the order parameter theory Sopl'l', ^n], has the 
McMillan-dcGennes form 

Sop = \N\^ J dvdt({dtVNf - Kiid^^f - K2(dy^N)^) 

+ J drdt(^\dt^\^-Cy\dy<P\^~C^\{d^~i^^N)\^-AcDwm^-ucDwm 

( 

where \N\ is the amplitude of the nematic order parameter, Ki and K2 are 
the two Franck constants (which were discussed before), Cx and Cy are the 
stiffnesses of the CDW order parameter along the x and y directions, Q is the 
modulus of the CDW ordering wavevector, Acdw and ucdw are parame- 
ters of the Landau theory that control the location of the CDW transition 
{Acdw = 0) and stability. Here we have assumed a stripe state, a unidirec- 
tional CDW, with its ordering wavevector along the x direction. We have also 
assumed z = 1 ( "relativistic" ) quantum dynamics which would be natural 
for an insulating system. 

The fermionic contribution has two parts. One part of the fermionic action, 
SflW\, where '0 is the quasiparticle Fermi field (we are omitting the spin in- 
dices), describes a conventional Fermi liquid, i.e. the quasiparticle spectrum 
with a Fermi surface of characteristic Fermi wavevector kp^ and the quasipar- 
ticle interactions given in terms of Landau parameters. What will matter to 
us is the coupling between the fermionic quasiparticles and the nematic order 
parameter (the complex director field N), and the CDW order parameter ^, 
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Sint^9N J drdt {Q2{r,t)N^{r,t) +h.c.)+gcDW J drdt {ncDw{r,t)<p1 {r,t) + h.c. 

(74) 

where qn and gcDW are two coupling constants and , as before, 

Q2{r,t)^^\r,t){d^+idyf4'{r,t) (75) 
is the nematic order parameter (in terms of quasiparticle Fermi fields), and 

ncDwiq,uj)^ J dkdn tp''{k + q + Cl,u} + n)tp{k,Q) (76) 

is the CDW order parameter (also in terms of the quasiparticle Fermi fields.) 

This theory has two phases: a) a nematic phase, for Acdw > 0, where 
)<?(= 0, and b) a CDW phase, for Acdw < 0, where {<P) ^= 0, separated 
by a QCP at Acdw = 0- In the nematic ("normal") phase the only low 
energy degrees of freedom are the (overdamped) fluctuations of the nematic 
Goldstone mode ifpf, and nematic susceptibility in the absence of lattice 
effects (which render it gapped otherwise) 

^^^^1' ^ glN{0) sin2(20,) - K,ql - K,ql) ^^^^ 

where sm{2(f>q) = 2qxqy/q^ and A^(0) is the quasiparticle density of states at 
the FS. 

We will consider here the simpler case in which the CDW ordering wavevec- 
tor obeys Q < 2kF (see the general discussion in Ref. |195| .) In this case one 
can see that the main effect of the coupling to the quasiparticlcs (aside from 
some finite renormalizations of parameters) is to change the dynamics of 
the CDW order parameter due to the effects of Landau damping. The total 
effective action in this case becomes 

/^CozM|<?(q,c.)|^ 

drdt (^Cy\dy^\^ + (^d^ ~ i^'fN^ + AcDwm^ + ucDwm 
+ (i^o^sin^(20,)-A-ig^-i^2<?^)|^^(q,a.)p 



(78) 



where Co ^ glow^ = 9l\N?N{Q) and Ki,2 = gl\N\^N{Q)K^,2. 

By inspecting Eg. ([75)) one sees that at Acdw = 0, as before the nematic 
Goldstone fluctuations have z — ?> (provided they remain gapless), and the 
CDW fluctuations have z = 2. Thus the nematic Goldstone modes dominate 
the dynamics at the nematic-CDW QCP. Even if the nematic Goldstone 
modes were to become gapped (by the lattice anisotropy), the QCP now will 
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have z = 1 (due to Landau damping) instead of z = 1 as in the "pure" order 
parameter theory. In both cases, the nematic Goldstone mode and the CDW 
order parameter fluctuations effectively decouple in the nematic phase. The 
result is that the nematic phase has relatively low energy CDW fluctuations 
with a dynamical susceptibility 

zCo|a;| - C^^qi - Cyq^ - Acdw 

(79) 

In other terms, as the QCP is approached, the nematic phase exhibits 
low energy CDW fluctuations that would show up in low energy inelastic 
scattering experiments much in the same was as the observed fluctuating 
stripes do in inelastic neutron scattering experiments in the high temper- 
ature superconductors . 1 1 1 7] As we saw before, a regime with "fluctuating" 
CDW (stripe) order is a nematic. 

A simple scaling analysis of the effective action of Eq. ((78| shows that, since 
2; > 1 at this QCP, the coupling between the nematic Goldstone mode tpN 
and the CDW order parameter <P is actually irrelevant. In contrast, in the 
(classical) case it is a marginally relevant perturbation leading to a fluctuation 
induced first order transition. [5H [53] Thus, this "generalized McMillan-de 
Gennes" theory has a continuous (quantum) phase transition which, possibly, 
may become weakly first order at finite temperature. 

This is not to say, however, that the stripe-nematic quantum phase transi- 
tion is necessarily continuous. In Ref. |I95| it is shown that the nature of the 
quantum phase transition depends on the relation between the ordering wave 
vector Q and the Fermi wave vector kp. For |Q| < 2fci? the transition is con- 
tinuous and has dynamical scaling z = 2. Instead, for |Q| = 2fci? it depends 
on whether |Q| is commensurate or incommensurate with the underlying 
lattice: for the incommensurate case the transition is (fluctuation induced) 
first order (consistent with the results of Rcf. [5]) but it is continuous for the 
commensurate case with z = 2 and anisotropic scaling in space. 

As in the case of the Pomeranchuk transition, the quasiparticles are effec- 
tively destroyed at the stripe-nematic QCP as well, indeed, already to order 
one loop it is found |195j that the quasiparticle scattering rate scales with 
frequency as S (cj) cx \/|ti7f, signaling a breakdown of Fermi liquid theory. 
As in our discussion of the nematic-FL QCP, this behavior must be taken as 
an indication of a breakdown of perturbation theory and not as the putative 
ultimate quantum critical behavior, which remains to be understood. 

In the quasiparticle picture we are using, the stripe state is similar to a 
CDW. Indeed, in the broken symmetry state the Fermi surface of the nematic 
is reconstructed leading to the formation of fermion pockets. As we noted 
above, we have not however assumed a rigid connection between the the 
ordering wave vector and the Fermi surface and, in this sense, this si not a 
weak coupling state. Aside from that, in the presence of lattice pinning of 
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the ncmatic Goldstonc mode, the asymptotic lowc energy properties of the 
stripe state are similar to those of a CDW (for details, see Ref. [195| .) 



9 Outlook 

In these lectures we have covered a wide range of material on the theory of 
electronic liquid crystal phases and on the experimental evidence for them. As 
it is clear these lectures have a particular viewpoint, developed during the past 
decade in close collaboration with Steven Kivelson. I have tried, primarily at 
the level of citations as well an on numerous caveats, to make it clear that 
there are many important unsolved and still poorly understood questions that 
(at present) allow for more than one answer. It is a problem that requires 
the development of many points of view which eventually complement each 
other. 

Several major problems remain open. One of them, in my view the most 
pressing one, is to establish the relation between the existence of these phases 
(stripes, nematics, etc.) and the mechanism(s) of high temperature supercon- 
ductivity. In my opinion there is mounting experimental evidence, some of 
which I discussed here, that strongly suggests the existence of a close and 
probably unavoidable connection. A question that deserves more considera- 
tion is the particular connection between nematic order and superconductiv- 
ity. Superficially these two issues would seem quite orthogonal to each other. 
Indeed, it is hard to see any connection within the context of a weak coupling 
theory. However if the nematic order arises from melting a stripe state which 
has a spin gap (such as the pair density wave state we discussed in these 
lectures) it is quite likely that a close connection may actually exist and be 
related. The current experimental evidence suggests such a relation. 

Another key theoretical question that is wide open is to develop a more 
microscopic theory of the pair density wave state. In spite with the formal 
analogy with the Larkin-Ovchinnikov state, it seems very unlikely that a a 
"straight BCS approach" would be successful in this case. This state seems 
to have a strong coupling character. 

As it must be apparent from the presentation of these lectures, the theory 
that has been done (and that is being done now) is for the most part quite 
phenomenological in character. There are very few "rigorous" results on the 
existence of these phases in strongly correlated systems. The notable excep- 
tion are the arguments we presented for the existence of nematic order in the 
strong coupling limit of the Emery model. Clearly more work is needed. 
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